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Abstract. The distribution of stellar masses that form together, the initial 
mass function (IMF), is one of the most important astrophysical distribution 
functions. The determination of the IMF is a very difficult problem because 
stellar masses cannot be measured directly and because observations usually 
cannot assess all stars in a population requiring elaborate bias corrections. Nev- 
ertheless, impressive advances have been achieved during the last decade, such 
that the shape of the IMF is reasonably well understood from low-mass brown 
dwarfs (BDs) to very massive stars. The case can be made for a rather universal 
form that can be well approximated by a two-part power-law function in the 
stellar regime. However, there exists a possible hint for a systematic variation 
with metallicity. From very elaborate observational surveys a picture is emerging 
according to which the binary properties of very-low-mass stars (VLMSs) and 
BDs may be fundamentally different from those of late-type stars implying the 
probable existence of a discontinuity in the IMF, but the surveys also appear to 
suggest the number of BDs per star to be independent of the physical conditions 
of current Galactic star formation. Star-burst clusters and thus globular cluster 
may, however, have a much larger abundance of BDs. Very recent advances have 
allowed the measurement of the physical upper stellar mass limit, which also ap- 
pears to be disconcertingly robust to variations in metallicity. Furthermore, it 
now appears that star clusters are formed in a rather organised fashion from low- 
to high stellar masses, such that the most-massive stars just forming terminate 
further star-formation within the particular cluster. Populations formed from 
many star clusters, composite populations, would then have steeper IMFs (fewer 
massive stars per low-mass star) than the simple populations in the constituent 
clusters. A near invariant star-cluster mass function implies the maximal cluster 
mass to correlate with the galaxy-wide star-formation rate. This then leads to 
the result that the composite-stellar IMFs vary in dependence of galaxy type, 
with potentially dramatic implications for theories of galaxy formation and evo- 
lution. 



1. Introduction 

The stellar initial mass function (IMF) is perhaps the most important macro- 
scopic distribution function in astrophysics, because it defines the mass-spectrum 
of stars born together. The other fundamental function is the star-formation 
history of a stellar system. They are connected through complex microscopic 
physical processes that defy detailed treatment on galactic scales. Together they 
contain the essential information on the transformation of dark gas to shining 
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stars and the spectral energy distribution thereof. They also contain the es- 
sential information on the cycle of matter, which fraction of it is locked up for 
"ever" in feeble stars and sub-stellar objects, and how much of it is returned en- 
riched with higher chemical elements to the interstellar medium or atmosphere 
of a galaxy. 

Given the importance of the IMF a huge research effort h as been invested 
to distill its shape and variability. The seminal contribution bv lSalpeter ( 19551 ) 



whilst staying in Canberra first described the IMF as a power-law, dN = 
£(m) dm = km~ a = dN, where dN is the number of stars in the mass in- 
terval m, m + dm and k is the normalisation constant. By modelling the spatial 
distribution of the then observed stars with assumptions on the star-formation 
rate, Galactic-disk structure and stellar evolution time-scales, Salpeter arrived 
at the power-law index a = 2.35 for OA^ui/Mq ^ 10, which today is known as 
the "Salpeter IMF". 

This IMF form implies a diverging mass density for m — > 0, which was 
interesting since dark matter was speculated, until the early 1990's, to possibly 
be made-up of faint stars or sub-stellar objects. Studies of the stellar velocities 
in the solar-neighbourhood also implied a large amoun t of missing, or dark, 
mass in the disk of the Milky Way (MW) ( Bahcalllll984l ). Careful compilation 



i n Heidelberg of the Gli ese Catalogue of Nearby Stars beginning in the 1960's 
( Jahreifi Sz Wielenl [19 971*1. and the application at the beginning of the 1980's 



of an innovative photographic pe ncil-beam survey-techn ique reaching deep into 
the Galactic field in Edinburgh bv lReid Gilmorg ( 1982I ) significantly improved 



knowledge of the space density of VLMSs (0.072 < m/M & < 0.14). 

Major studies extending Salpeter's work to lower and larger masses followed, 
showing that the mass function (MF) of Galactic-field stars turns over below one 
solar mass thus avoiding the divergence. Since stars with masses m^0.8M Q 
do not evolve significantly over the age of the Galact ic disk , the MF equals 
the IMF for these. While the work of Miller Scald ( 19791) re li ed on using 



the nearby stellar sample to define the IMF for m < 1 Mm. IScalol (Il986l ) relied 
mostly on a more recent deep pencil-beam star-count survey. IScalol (|1986l ) stands 
out as the most thorough and comprehensive analysis of the IMF in existence, 
laying down notation and ideas still in use today. The form of the IMF for low- 
mass stars was revised in the early 1990's in Cambridge, especially through the 
quantification of significant non-linearities in the stellar mass -luminosity relation 
and evaluat ion of the bias due to unresolved binary systems ( Kroupa et al.|[l990l . 



119911 . 119931 ). On the one hand this work lead to a detailed understanding of 
the shape of the stellar luminosity function (LF) in terms of stellar physics, 
and o n the other t he di fference between the results obtained by iMiller Sz Scald 
(|1979T ) and IScald (jl986h for low -mass stars was resolved by this work through 



rigorous modelling of all biases affecting local trigonometric-based and distant 
photometric-parallax-based surveys, such as come from an intrinsic metallicity 
scatter, evolution along the main sequence and contraction to the main sequence. 
In doing so this work also included an up dated local st e llar s ample and the then 
best-available deep pencil-beam survey bv lStobie et al.l ( 19891 ). As such it stands 



1 The latest version of the catalogue can be found at http://www.ari.uni-heidelberg.de/aricns/ 



while http://www.nstars.nau.edu/ contains the Nearby Stars (NStars) database 
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unique today as being the only rigorous analysis of the late-type-star MF using 
simultaneously both the nearby trigonometric parallax and the far, pencil-beam 
star-count data. This study was further ext en ded to an ana l ysis o f other ground- 
based pencil-beam surveys (lKroupal[l995al ). IZheng et all pOOll ) employed the 
HST to measure the LF through the entire thickness of the MW disk finding 
excellent agreement wit h the pencil-b e am (p hotometric-parallax) LFs (Fig. [U] 
below). The results of iKroupa etail (jl993T ) (a rj 1.3,0.08 - 0.5 M ), were 
confirmed by iReid et al.1 (|2002l ) using updated local star-counts that included 
Hipparcos parallax data. Indeed, the continued long-term observational effort 
on nearby stars by Neill Reid, John Gizis and collaborators forms one of the 
very major pillars of modern IMF work; continued discussion of controversial 
interpretations have much improved and sharpened our general understanding 
of the issues. A re-analysis of the nearby mass function of stars in terms of 
a log-normal form in the mass range 0.07 — 1 M & was provided by Chabrier, 
finding agreement to the deep HST star-count data once unresolved multiple 
stars and a metal -deficient colour- magnitude relation for thick-disk M dwarfs 
are accounted for ( Chabrienl2003bl ). 

The Galactic-field IMF for 0.08 < vi/Mq < 1 can thus be regarded as being 
reasonably well-constrained, but some unresolved issues nevertheless remain. 
The exact form of the IMF is still under dispute. The mathematical necessity for 
the IMF to turn-o v er at some small mass and the corresponding empirical result 
bv lMiller Sz Scalo ( 19791 ) lead to the development of a theory of th e IMF based 
on gravitational in stabilities in the turbulent inter-stellar me dium (lFlecklll982l ; 
iFerrini et al. 19831 ) and random hierarchical fragmentation (Zi nneckerl I'l984). 
thereby accounting for the approximately log-normal shape of the IMF. While a 
log-normal form has often been adopted given th e turnover near a few tenths of a 
solar mass dMiller fc Scaldll979l:IChabrier{ l2001) and some theoretical ideas (e.g. 
Adams &: Fatuzzolll996l ; IChabrierl l2003al ) . it has never been demonstrated that 



the log-normal form is consistent with the sharp peak in the stellar luminosity 
function (LF) near absolute visual and I-band magnitudes My ~ 12, Mi ~ 8.5, 
respectively. A power-law description together with a semi-empirical mass- 
lumi nosity relation has, in contrast, been found to fit the nearby and deep LFs 
well ( Kroupa et al.|[l993l ; lKroupalll995al ). Furthermore, most approaches have 
relied on using parametrised formulae, such as the multi-power-law form or the 
log-normal form, and so it would be important to develop a non-parametric 
method that allows inclusion of the biases through unresolved multiple systems. 
Such a method would be very useful for studying the shape of the IMF in star 
clusters. 

In contrast to the Galactic-field sample, where stars of many ages and 
metallicities are mixed, star clusters offer the advantage that the stars have 
the same age and metallicity and distance. And so a very large effort has 
been invested t o try to extract the IMF from open (ISanner Geffertl l200ll 
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Piskunov et al. 20041 ; Moraux et al.l 12004 e.g.) and embedded (jHillenbran 
19971: iMuench et al.ll2002l : lLuhmanll2004l. e.g.) clusters, as well as associations 
(|Gouliermis et al.l I2005T T 3ere the continued methodologically consistent ob- 



servational work on a number of different very young populations by Kevin 
Luhman and collaborators has had a major impact on our understanding of 
the IMF at low masses. On the theoretical side, the ever-improving modelling 
of stellar and BD atmospheres being pushed forward with excellent results no- 
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tably by the Lyon group (Isabelle Baraffe and Gilles Chabrier), has allowed 
consistently better constraints on the faint-star MF by a wide variety of ob- 
servational surveys Q. Furthermore, the development of high-precision iV-body 
codes that rely on complex mathematical and a lgorithmic pro cedures through 
the work of Sverre Aarseth and Seppo Mikkola ( Aarseth 19991 ) and others has 
led to important progress on understanding the variation of the dynamical 
propertie s of st ellar populations in i ndividual clusters and the Galactic field 
(lKroupalll995dl: Ide La Fuente Marcos! 1 19971: lKroupall2001bl: iKroupa et all 120011: 
Portegies Zwart et al.ll200ll , l2002l ; lBaumgardt fc Makindl2003l ; |Moraux et alj|200 



104; 



Hurley et al.ll2005l ). In general, the MF found for clusters is consistent with the 
Galactic-field IMF for m < IMq, but some problems remain. While a huge 
progress has been achieved in modelling cluster populations over time, it is for 
exampl e still unclear as to wh y open clusters have a significant deficit of white 
dwarfs dFellhauer et al.l 12003). A further pr oblem is posed by the rapid and 



violent early dynam ical evolution of clusters ( Lada et al. Il984 ; iGoodwinl 1997 ; 
Kroupa et aD 120011 : iGever & Burkertl 120011 : iMoraux et al.ll2004l : iKroupal 120051 ) 
and the associated loss of a large fraction of the cluster population, and due to 
the density-dependent disruption of primordial binary systems. Young clusters 
have thu s undergone a highly complex dynamical evolution which continues into 
old age ( Baumgardt &: Makind 12003) and are therefore subject to biases that 
can only be studied effectively with full-scale iV-body methods, thus imposing 
a complexity of analysis tha t surpasses th at for the Galactic-field sample. 
For massive stars, the IScald (1986) 

IMF is based on a combination of 
Galactic-field star-counts and OB association data and has a slope a ~ 2.7 
for m£ 2M(?, with mu c h unc ertainty for m<;10 Mq. The previous determina- 
tion bv lMiller fc Scaio ( 19791 ) also argued for a relatively steep field-IMF with 
a = 2.5,1 < m/MQ < 10 and a = 3.3, m > 10 Mq. However, Massey's work 
at Tucson demonstrated, t hrough extens ive spectroscopic classification, that 
a = 2.35 ± 0.1,m;^ 10 M ( Massevl Il998l ). for a large variety of physical en- 
vironments, namely OB associations and dense clusters for populations with 
metallicity ranging from near-solar-abundance to about l/10th metal abun- 
dance (Fig. [T]). Consequently, above about 0.5 Mq the empirical IMF can 
be described well by a single power-law form with the Salpeter/Massey index, 
a = 2. 35, therewith being remarkably invariant 



While the case is often ma de that star 
abundant in massive stars (e.g. Eisenhauer 



tion mimics just such an effect (jBoilv et al 



aurst s may have an IMF that is over- 
time-dependent mass segrega- 
20051 ) . Local well-resolved star-burst 



clusters also do not support top-heavy IMFs. Notable well-studied cases with 
masses in the range 1O 4_5 M0 and central densities near 10 5 stars/pc 3 are the 



2 Here it should be emphasised and acknowledged that the intensive and highly fruitful discourse 
between Guenther Wuchterl and the Lyon group has led to the important understanding that 
the classica l evolution tracks computed in Lyon and by others are unreliable for ages less than 
a few Myr (| Wuchterl fc Tscharnuter]|2003D . 

3 Note that!Scalo ( 1998) emphasises that the IMF remains poorly constrained owing to the small 
number of massive stars in any one sample. This is a true albeit conservative stand-point, and 
the present author prefers to accept Massey's result as a working hypothesis. This hypothesis 
that there exists a universal parent distribution func tion is tested in individual clusters and 
OB associations for possible significant deviations (cf. iPflamm- Altenburg fc Kroupd 120061 ). 
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30 Dor cluster (R136) in the LMC, NGC 3603 in the MW, and the Arches cluster 
near the Galactic centre. The 30 Dor star-burst cluster (NGC 2070) has been 
found bv lMassev Hunted (119981) to have a Salpeter MF for 2.8 <m/M e < 120, 
which is confirmed by Selman et al. ( 19991 ) who apply corrections for variable 
reddening. But they also note that the c ore of NGC 2070 is clearly mass- 
segregated with a flatter MF. For NGC 3603. ISung & Besselll (|2004l ) find a highly 
mass-segregated system with a flat MF in the core, but the MF beyond the core 
has a = 1.9 ±0.1, m ^ 1.6 M@, thus again being consistent with a Salpeter value. 
The Arches cluster is situated at a projected distance of 25 pc from the Galactic 
centre and i s therefore in a ver y ex otic environment be ing strongly influenced by 
tidal forces. iFiger et all ( 19991 ) and lStolte et al.l ((2002) fin d the innermost r egion 
of the cluster to have a ss 1, but at larger radii a « 2.7. iKim et all ( 20001 ) and 



Portedes Zwart et al.1 (|2002h study the evolution of the Arches using iV-body 



modelling and find that the cluster evolves rapidly, loosing memory of its birth 
configuration within about 1 Myr. In particular, they find that the IMF in the 
Arches may have been quite normal. 

This short discussion thus indicates that the evidence for a top-heavy IMF 
is not strong in well-resolved star-burst clusters, and that dynamical evolution 
of the clusters needs to be modelled in detail to under stand p o ssibly devia nt 
observed IMF's. A highly interesting finding reported by iMassevl (1998, 20031) is 
that the OB stellar population found in the field of the LMC has a very steep 
MF, a ~ 4.5, which can be interpreted to be due to the preferred formation 
of small groups or even isolated O and B-type stars. Another interpetation 
that would not need to resort to exotic star-fomration events is the dynamical 
ejection of OB stars from dynamically unstable cores of young clusters. This 
may lead to such a steep IMF because simple esti mates show that prefere ntlially 
the less-massive members of a core are ejected (IClarke fc Pringldll992D . This 
process needs to be studied using fully-consistent and high-precission N— body 
modelling of youing clusters, to see if the observed distribution of field-OB stars 
can be accounted for with this process alone, or if indeed an exotic star-formation 
mode needs to be invoked to explain some of the observations. 

The "secret" ingredient for uncovering the true nature of an IMF is thus the 
high-precision N-body experiment , which therewith achieves an importance well 
beyon d merely the methodological \Kim et ~al\\200d ; \Pflamm-Altenbura & Kroupd 
MM, e.g.). ' 

Although we have essentially no knowledge about the low-mass end of the 
IMF in star-burst clusters, due to their rarity and thus generally large dis- 
tances, today we do know that the MF of low-mass stars in globular clusters, 
that presumably formed as star-burst clusters, is quite similar to that of the 
solar- neighbourhood and in youn g open clusters. This sugge sts a remarkable 
invar iance of the low- mass IMF ( de Marchi Sz Paresce 1995al lbl; iParesce et al .1 
19951 . e.g.). The case has been made for a systematic variation with metallicity 



in the sense that the metal-poorer and older pop ulations may h ave flatter MFs 
as expected from simple Jeans-mass arguments (lKroupd[2TOlbl . eq. [22] below) , 
but this suggestion remains speculative r ather than conclusive because stellar- 
dyna mical processes skew observed MFs ( Kroupa 2001bl ; Baumgardt &: Makinol 
20031 ) . 
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The MF for BDs is even shallower, as two recent studies of the available 
solar-neighbourhood star-count data demons trate, these constraints being con - 
sistent with those from young star cluste rs (lChabrierll2002l; lAllen et al.l 12 005). 



From their non-s tellar binary properties (Martin et al. 20031 ; Bouv et al. 20031 : 



Close et alJl2003h it is emerging, however, that BDs may not be a continuation 
of the stellar population. Although their formation is linked to that of their 
sibling stars, the IMF is likely to be discontinuous near the hydrogen-burning 
mass limit , and the overall result a ppears to be that about one BD is born per 
five stars ( Kroupa Sz Bouvierl 120031 . § 15.2.1 below). We have thus had to realise 
that BDs contribute a negligible dynamical mass to any stellar population. Faint 
main-sequence stars account for most of the mass in the MW disk, and the need 
for dark matter in the MW disk also disappeared as improved kinematica l data 
and analysis thereof became available (|Kuiikedll99l : lFlvnn fc Fuchsll 19941 1. Fur- 
thermore, the suggestion has been made by two independent teams, based on 
an extremely high-quality spectral analysis of solar-neighbourhood stars and by 
accounting for the relative death rate of stars, that the G alactic thick disk may 
be substantially more massive than is currently thought ( Soubiran et al.ll2003l ; 
Fuhrmann 20041 ) . This would further reduce the mass discrepancy evident in 



the solar motion about the Galaxy. 

Five very recent further devel opments have enrich e d the study of the mass- 
distribution of stars appreciably. IWeidner Sz Kroupa ( 20041 ) demonstrated for 
the first time that the absence of stars more massive than about 150 Mq in R136 
in the LMC must imply a physical mass limit to stars near 150 Mq, with high 
statistical significance, unless a ^2.8 for m^lM . Such a steep IMF may be 
the true mass-distribution if unresolved multiple systems are corrected for, but 
an affirmation of this and thus the negati on of the ph ysical stellar mass limit 
awaits a detailed investigation of this issue. iFigerl ( 20051 ) followed by performing 
the same analysis with hi s HST data on t he Ar ches cluster confirming the same 
physical mass limit, and lOev Sz Clarke! ( 20051 ) published a statistical analysis 
of a number of OB associations and clusters again confirming that a physical 
mass limit exists in th e mass range 12 — 200 Mq. Thus, stellar masses are 
limited near 0.072 M Q (jChabrieri 120033 ) and near 150 Mq. The other notion 
recently re-affirmed is that star clusters appear to limit the masses of the most 
massive stars in them in such a wa y that more massive clust ers permit higher 
masses of their most massive stars ( Weidner &: Kroupa 2006). The third very 
recent insight, resulting from the above two, is that the composite IMFs of 
whole galaxies must be steeper than the stellar IMF for m^l Mq, because 
composite populations result from the addition of the stellar IMFs in the star 
clusters which themselves are distr i buted according to a powe r-law cluster mass 



function ( Kroupa k, Weidnei 20031 ; Weidner Kroupall2005l ). This notion had 

3n| i 



be en raised already bvlVanbevererl t 19821 ) and a similar ansatz has been followed 
by lOev &; Clarkd ( 1998! ) for the construction of HII LFs in galaxies, and the 
recent advances on this issue thus suggest that the discrepancy between Scalo's 
field-star IM F and Massey's stellar IM F in star clusters may have found a natural 
resol ution ( Kroupa Sz Weidner 20031 ). The fourth recent insight ( Weidner et al 



2004) is that the star-formation rate of a galaxy dictates the mass of the most 
massive star cluster formi ng within it. This then implies the fifth recent result 
( Weidner Sz Krou pa 2006) that galaxies have different integrated (composite) 
IMFs. According to this picture, galaxies with a small stellar mass ( ^ 10 Mq) 
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may, over their life-times, be significantly deficient in stars more massive than a 
few M Q . 

Given these results, it thus becomes necessary to distinguish between simple 
stellar populations and composite populations. A simple population is found in 
a star cluster and consists of stars of equal age and metallicity. A composite 
population consists of more than one cluster. The stellar IMF refers to the 
IMF of stars in a simple population, while the composite IMF, or the integrated 
galactic IMF (IGIMF) is the IMF of a composite population, i.e. a population 
composed of many star clusters, most of which may be dissolved. 

This treatise attempts to provide an overview of the general methods used 
to derive the IMF with special attention on the pitfalls that are typically en- 
countered. The somewhat uncomfortable result until now is that there is no 
reliable, or confirmed, evidence for stellar-IMF differences in simple populations 
with different physical properties. While this must be recognised as being a 
problem, because elementary arguments would imply variations with different 
conditions, it does ease modelling of simple stellar systems. While the IMF for 
0.08 m/M & ^ 1 M© is by now quite well constrained, for massive stars the ef- 
fects of unresolved multiple systems remains unknown^. Corrections for these 
may significantly affect the shape of the IMF with rather major implications for 
much of extragalactic and Galactic astrophysics. Understanding these correc- 
tions is thus of major importance. The very recent realisation that star-clusters 
limit the mass spectrum of their stars has interesting implications for the for- 
mation of stars in a cluster and leads to the insight that composite populations 
must show IMFs that differ from the stellar IMF in each cluster. A detailed 
analysis indicates that the composite IMF, i.e. the IGIMF, should vary between 

galaxies. With this finale, this treatise reaches the cosm ologica l arena . 

Ot her relatively rec ent reviews of the IMF are by IScald (|1998T ). iKroupal 

( 20021 ). IChabrieil ( 2003al ). The proceedings of the "38th Herstmonc eux Con- 
ference on the Stellar Initial Mass Function" ( Gilmore &; Howell Il998l ) and the 
proc eedings of the "IMF5 0" meeting in celebration of Ed Salpeter's 80th birth- 
day (ICorbelli et al.1 12005T ) contain a wealth of important contributions to the 
field. 



2. Some essentials 

Assuming all binary and higher-order stellar systems can be resolved into indi- 
vidual stars in some population such as the solar neighbourhood and that only 
main-sequence stars are selected for, then the number of single stars per pc 3 in 
the mass interval m to m + dm is dN = dm, where H(m) is the present- 

day mass function (PDMF). The number of single stars per pc 3 in the absolute 
V-band magnitude interval M P to M P + dM P is dN = -^>{M P )dM P , where 
$(Mp) is the stellar luminosity function (LF) which is constructed by counting 
the number of stars in the survey volume per magnitude interval, and P signifies 



4 see also footnote [3] 
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an observational photometric pass-band such as the V- or /-band. Thus 

Note that the the minus sign comes-in because increasing mass leads to decreas- 
ing magnitudes, and that the LF constructed in one photometric pass band P 
can be transformed into another band P' by 

. . dN dM p/ T/ „, s dM P : 

if the function Mpi = fn(Mp) is known. 

Since the derivative of the stellar mass-luminosity relation (MLR), m{Mp) = 
m(Mp, Z,t,s), enters the calculation of the MF, any uncertainties in stellar 
structure and evolution theory on the one hand (if a theoretical MLR is relied 
upon) or in observational ML-data on the other hand, will be magnified accord- 
ingly. This problem cannot be avoided if the mass function is constructed by 
converting the observed stellar luminosities one-by-one to stellar masses using 
the MLR and then binning the masses, because the derivative of the MLR never- 
theless creeps-in through the binning process, because equal luminosity intervals 
are not mapped into equal mass intervals. The dependence of the MLR on the 
star's chemical composition, Z, it's age, r, and spin vector s, is explicitly stated 
here, since stars with fewer metals than the Sun are brighter (lower opacity), 
main-sequence stars brighten with time and loose mass, and rotating stars are 
dimmer because of the reduced internal pressure. Mass loss and rotation are 



signi ficant factors for intermediate and especially high-mass stars (jPennv et al 
12001] ). 

The IMF follows by correcting the observed number of main sequence stars 
for the number of stars that have evolved off the main sequence. Defining t = 
to be the time when the Galaxy that now has an age t = tq began forming, the 
number of stars per pc 3 in the mass interval m, m + dm that form in the time 
interval t, t+dt is dN = £(m, t) dm x b{t) dt, where the expected time-dependence 
of the IMF is explicitly stated (§ [8]), and where b(t) is the normalised SFH, 
(1/tq) Jq TG b(t) dt = 1. Stars that have main-sequence life-times r(m) < tq leave 
the stellar population unless they were born during the most recent time interval 
r(m). The number density of such stars with masses in the range m,m + dm 
still on the main sequence and the total number density of stars with r(m) > tq, 
are, respectively, 

. . „. , 1 f r° , Mt)dt , rim) < tq, , , 

Z(m) = £(m)— x < 7£T 7 T (™K , ( (3) 

V ' ^ t G \ £ G b(t)dt , T{m) >r G , V ; 

where the time-averaged IMF, £(m), has now been defined. Thus, for low-mass 
stars E = £, while for a sub-population of massive stars that has an age At <^ tq, 
E = £ (At/rc) for those stars of mass m for which the main-sequence life-time 
r(m) > At, indicating how an observed high-mass IMF in an OB association, 
for example, has to be scaled to the Galactic-field IMF for low-mass stars0 In 



5 assuming continuity of the IMF 
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this case the different spatial distribution via different disk-scale heights of old 
and young stars also needs to be taken into account, w hich is done globally 
by calculat ing the stellar surface density in the MW disk (jMiller Scalolll979i ; 
Scald Il98d ). Thus we can see that joining the cumulative low-mass star counts 



to the snap-shot view of the massive-star IMF is non-trivial and affects the 
shape of the IMF in the notorious mass range ~ 0.8 — 3M , where the main- 
sequence life-times are comparable to the age of the MW disk. In a star cluster 
or association with an age r c i <C tq, t c \ replaces tq in eq. [3l Examples of the 
time-modulation of the IMF are b{t) = 1 (constant star-formation rate) or a 
Dirac-delta function, b{t) = t c \ x S(t — to) (all stars formed at the same time to). 
Often used is the "logarithmic mass function" (Table below), 

£ L (m) = (mlnlO) £(m), (4) 

where dN = £,h( m ) dim is the number of stars with mass in the interval Im, lm + 
dim {lm = log lo m)0. 



3. The massive stars 

Studying the distribution of massive stars is complicated because they radiate 
most of their energy at far-UV wavelengths that are not accessible from Earth, 
and through their short main-sequence life-times, r, that remove them from 
star-count surveys. For example, a 8 5 star cannot b e distinguished from a 
40 M & star on the basis of My alone (|MassevHl998l . 120031 ') . Constructing *(My) 
in order to arrive at E(m) for a mixed-age population does not lead to success 
if optical or even UV-bands are used. Instead, spectral classification and broad- 
band photometry for estimation of the reddening on a star-by-star basis has 
to be performed to measure the effective temperature, T e fj, and the bolometric 
magnitude, Mbd, from which m is obtained allowing the construction of S(m) 
directly (whereby $(Mb i) and H(m) are related by eq. [I]). Having obtained 
H(m) for a population under study, the I MF follows by applying eq. [31 A more 
straight-forward method often used (e.g. lHillenbrandlll997l ; lMassevill998l . 120031 ) 



is to evolve each measured stellar mass to its initial value using theoretical stellar 
evol ution tracks, a nd to construct £(to) from this set of masses. 

Im assevl (|l003) stresses that studies that only rely on broad-band optical 



photometry consistently arrive at IMFs that are significantly steeper with 03 ~ 3 
(see eq. [201 below), rather than 03 = 2.2 ± 0.1 consistently found for a wide 
range of stellar populations. Indeed, the application of the same methodology 
by Massey on a number of young populations of different metallicity and density 
shows a remarkable uniformity of the IMF above about 10 M@ (Fig. [1]). 

The available IMF measurements do not take into account the bias through 
unresolved systems which may be substantial since the proportion of multi- 
ple stars is higher fo r massive stars than for low-mass Galactic-field stars (e.g. 
Duchene et al.ll200ll ). For example, in the Orion Nebula Cluster (ONC) each 



massive star has, o n average, 1.5 compan i ons dPreibisch et~aL 19991 ). while in the 



cluster NGC 6231 iGarcia fc Mermilliod! (|200l|) find that 80 % of all O stars are 



'Note that IScalol (| 19861 ) calls t;h(rn) the mass function and £(m) the moss spectrum. 
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Fi gure 1. The IMF slope T — I — a determined in a homogeneous manner 
by IMasseyl (|2003l ) for OB associations and clusters in the MW, LMC and 
SMC. The Small Magellanic Cloud (SMC) has a metallicity Z = 0.005, the 
Large Magellanic Cloud (LMC) has Z = 0.008 and the Milky Way (MW) has 
Z = 0.018 wi t hin a distance of 3 kpc around the Sun. With kind permission 
from IMasseyl (|2003l ). 



radial- velocity binaries. Sagar Sz Richtler ( 199ll ) estimate the binary-star bias 
on the IMF for stars in the mass range 2 to 14 M© assuming each star has one 
companion. The IMF would steepen to the Scalo value a ~ 2.7 for a measured 
a = 2.3. This correction, howev er, depends on the distribution of compan ions 
that is not yet known very well (|Preibisch et al.lll999l : [Duchene et al.ll200ll L A 
larger value, a 3 ± 0.1, is also suggested from a completely independent but 
indirect app roach relying on the distribution of ultr a-compact HII regions in 
the Galaxy (|Casassus et al.ll2000l : lOev fc ClarkdfTTwsL but this may be a result 
of the composite nature of galaxy- wide populations (§ \8$. Clearly, the effect 
multiplicity has on the massive-star IMF needs further exploration. 

Massive main-sequence stars have su bstantial winds flowing outwards with 
velocities of a few 100 to a few 1000 km/s ( Kudritzki Sz Pulsll200Cll ) . For example, 
10~ 6 5 < M < 10~ 6 M /yr for m = 35 M Q with r = 4.5 Myr (|Garcia-Segura et alj 
19961). and 10~ 5 6 < M < lO _5 - 8 M /yr for m = 60 M with r = 3.45 Myr 



( Garcia-Segura et al. 19961 ). More problematical is that stars form rapidly ro- 
tating and are sub-luminous as a result of reduced internal pressure. But they de- 
celerate significantly during their main-sequence life-time owing to the angular- 
momentum loss throug h their winds and becom e more luminous more rapidly 
than non-rotating stars ( Maeder &: Mevnet 2000l ). A comparison of such models 
is available in Figs. [2] and [3l Evidently, for ages less than 2.5 Myr the models 
deviate only by 5-13 % from each other in mass, luminosity or temperature 
( Weidner &; Kroupa 2006). Large deviations are evident for advanced stages of 
evolution though. 

The mass-luminosity relation for a population of stars that have a range 
of ages is therefore broadened making mass estimates from Mboi uncertain by 
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Figure 2. The mass-evolut ion of massive stars a ccording; to theoretical work: 
Solid lin es: Geneva models JSchaller et alJl!992f ). Dashed lines: Geneva with 
rotation (jMevnet fc Maederll2003l . (thick: none, medium: 300 km/s, thin: 500 
km/s;) Dotted lines: stellar evolution pa ckage fr om Hurley based on the P ols 
non-rotating models (|Hurlev et al.l f2000). From Wci dner fc Kroupal (|2006| ). 
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Figure 3. The luminosity-evolution of massive stars according to the- 
oretical work. Line-typ es and initial masses are as in Fig. [21 From 
IWeidner fc Kroupal (|2006l ) 
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up to 50 per cent ( Penny et al.ll200ll ). a bias that probably needs to be taken 
into account more thoroughly in the derivations of the IMF. Another problem 
is that m ^ 40 Mq stars may finish their assembly after burning a significant 
proportion of their c entral H so that a zero-a ge-main sequence may not ex- 
ist for massive stars ( Maeder &: Behrendl I20021 ). However, the agreement be- 
tween slowly-rotating tidally-locked massive O-t ype binaries with s tandard non- 
rotating theoretical stellar models is very good (|Pennv et al.ll200lh . 



3.1. The maximum stellar mass and cluster formation 



A brief history: The empirically determined range of stellar masses poses im- 
portant constraints on the physics of stellar formation, structure and stellar 
evolution, as well as on the feedback energy injected into a galaxy's atmosphere 
by a population of brand-new stars. The physical limit at low masses is now well 
established (§I5T|). and an upper mass limit appears to have been found recently. 

A theoreti cal physica l limit ation to stellar masses has been known since 
many decades. lEddington ( 19261 ) calculated the limit which is required to bal- 
ance radiation pressure and gravity, the Eddington limit: Le6&/Lq « 3.5 x 
10 4 th/Mq. Hydrostatic equilibrium will fail if a star of a certain mass m has a 
theoretical luminosity that exceeds this limit, which is the case for m^60 Mq. 
It is not clear if stars above this limit cannot exist, as massive stars are not 
fully radiative but have convective cores. But more massive stars will l oose 
material rapidly due to strong stellar winds. Schwarzschild &: Harm ( 19591 ) in- 
ferred a limit of « 60 M beyond which stars should be destro yed due to pulsa- 
tions . But later studie s suggested that these may be damped (jBeech fe Mitalas 
19941 ). IStothersI (|l992l ) showed that the limit increases to m max * « 120- 150 M 
for more recent Rogers-Iglesia opacities and for metallicities [Fe/H]~ 0. For 
[Fe/H]« -1, 

nimax* ~ 90 Mq . A larger physical mass limit at higher metallicity 
comes about because the stellar core is more compact, the pulsations driven by 
the core having a smaller amplitude, and because the opacities near the stellar 
boundary can change by larger factors than for more metal-poor stars during 
the heating and cooling phases of the pulsations thus damping the oscillations. 
Larger physical mass limits are thus allowed to reach pulsational instability. 

Related to the pulsational instability limit is the problem that radiation 
pressure also opposes accretion for proto-stars that are shining above the Ed- 
dington luminosity. Therefore the question remains how stars more massive 
than 60 Mq may be formed. Stellar formation models lead to a mass limit 
near 40 — 100 imposed by feedback on a spherical accre tion envelope (|Kahn 



19741 ; IWolfire fc Cassinelli|[l986l ; IWolfire k Cassinellil fl987l ) . Some observations 
suggest that star s may be accreting material in discs and not in spheres (e.g. 



Chini et al 



2004j). The higher density of the disc- material may be able to over- 
come the radiation at the equator of the proto-star. But it is unclear if the 
accretion-rate can be boosted above the mass-loss rate from stellar winds by 
this mechanism. Theoretical work on the formation of massive stars through 
disk-accretion with hig h accretion rates thereby allowing th ermal radiation to 
escape pole- wards (e.g. lNakanolll989l ; iJiiina Adams! Il996l ) indeed lessen the 
problem and allow stars with larger masses to form. 

Another solution proposed is the merging scenario. In this case massive 
stars form through the merging of intermediate-mass proto-stars in the cores 
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of dense stellar clusters driven by core-contraction due to very rapid accretion 
of gas with low specific angul ar momentum , thus again avoiding the theoretical 
feedback-induced mass limit (|Bonnell et al.lll998l : IStahler et al.) 120001 ) . It is un- 
clear though if the very large central densities required for this process to act 
are achieved in reality, but it should be kept in mind that an observable cluster 
is, by necessity, exposed from its natal cloud an d is therefore likely to be always 
observed in an expanding phase ( Kroupall2005l ). 

The search for a possible maximal stellar mass can only be performed in 
massive, star-burst clusters that contain sufficiently many stars to sample the 
stellar initial mass function beyond 100 Mq. Observationally, the existence of 
a finite physical stellar mass limit was not evident until very recently. Indeed, 
observations in the 1980's of R136 in the Large Magellanic Cloud (LMC) sug- 
gested this objec t to b e one single star with a mass of about 2000 — 3000 M . 
Weigelt fc Baier ( 19851 ) for the first time resolved the object into eight com- 
ponents using digital speckle interferometry, therewith proving that R136 is a 
massive star cluster rather than one single super-massive star . The evidence fo r 
any physical upper mass limit became very uncertain, and Elmegreen ( 19971 ) 
stated that "observational data on an upper mass cutoff are scarce, and it is 
not included in our model s [of the IMF from rando m sampling in a turbulent 
fractal cloud]". Although iMassev &: Hunter (119981 ) found stars in R136 with 
masses ranging up to 140 — 155 Mq. IMassev (|2003l ) explains that the observed 
limitation is statistical rather than physical. We refer to this a s the Massey 
assertion, i.e. that m max * = oo. Meanwhile, ISelman et al.l ( 19991 ) found, from 



their observations, a probable upper mass limit in the LMC near about 13 M^ 
but th ey did not evaluate the statistical significance of this suggestion. iFigerl 
(2003) discussed the apparent cut-off of the stellar mass-spectrum near 150 Mq 
in the Arches cluster near the Galactic centre, but again did not atta ch a statis- 
tical analysis of the significance of this observation. lElmegreenl (|2000l ) also noted 
that random sampling from an unlimited IMF for all star-forming regions in the 
Milky Way (MW) would lead to the prediction of stars with masses <L 1000 Mq, 
unless there is a rapid turn-down in the IMF beyond several hundred Mq. How- 
ever, he also stated that no upper mass limit to s tar formation has ever been 
observed, a view also emphasised by iLarsonl ( 20031 ). 

Thus, while theory clearly expected a physical stellar upper mass limit, 
the observational evidence in support of this was very unclear. This, however, 
changed in 2004. 



Empirical res ults Given the observed r ather sharp drop-off of the IMF in R136 
near 150 Mq, Weidner fe Kroupa ( 20041 ) studied the Massey assertion in some 

detail. 

R136 has an age < 2.5 Myr ( Massey Hunter 1998! ) which is young enough 
such that stellar evolution will not have removed stars through super nova ex- 
plosions. It has a metallicity of [Fe/H]sa — 0.5 dex (de Boer et al.|["l985l ). 

From the radial surface density profile ISelman et~ai1 ( 19991 ) estimated there 
to be 1350 stars with masses between 10 and 40 Mq within 20 pc of the 30 Do - 
rad us region, with i n the centre of which lies R136. IMassev Hunter] ( 19981 ) 
and ISelman et all (11999T ) found that the IMF can be well-approximated by a 
Salpeter power-law with exponent a = 2.35 for stars in the mass range 3 to 
120 M . This corresponds to 8000 stars with a total mass of 0.68 x 10 5 M 



0- 
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Extrapolating down to 0.1 M@ the cluster would contain 8 x 10 5 stars with a 
total mass of 2.8 x 10 5 M & . Using a standard IMF with a slope of a = 1.3 (in- 
stead of the Salpeter value of 2.35) between 0.1 and O.5M this would change 
to 3.4 x 10 5 stars with a combined mass of 2 x 10 5 M o , for an average mass 
of 0.61 M(D oyer the mass range 0.1 — 120_Mq . Based on the observations by 
Selman et aD ([1999ft . IWeidner &; Kroupal (l2004h assumed that R136 has a mass 
in the range 5 x 10 4 < M-rxzq/Mq < 2.5 x 10 5 . This mass range can be used to 
investigate the expected number of stars above mass m, 



N(> m) 



£(m ) dm' , 



(5) 



with the mass in stars of the whole (originally embedded) cluster being 



M ec i 



m £(m) dm' , 



(6) 



where m\ ow = 0.01 M & and m max * = oo (the Massey assertion). Here the as- 
sumption is that the cluster is still compact despite having-blown out its residual 
gas. There are two unknowns (N(> m) and k) that can be solved for using the 
two equations above. 

Using the standard stellar IMF (eq. 1201 below). N(> m) is plotted in Fig. H] 
for the two mass estimates of the cluster. The solid vert ical line indicates 
150 M ^, the approximate maximum mass observed in R136 (jMassev fc Hunter 
19981 ) . IWeidner fc Kroupal (120041 1 found that N(> 150 M 



i , 

ing if M ecl = 2.5 
M Prl = 5 x 10 4 M, 



40 stars are miss- 

: 10 5 M Q , while N(> 150 M ) = 10 stars are missing if 
q. The probability that no stars ar e observed although 10 are 



expected, assuming m max * = oo, is P = 4.5 x 10 -5 . IWeidner &; Kroupal ( 20041 ) 
concluded that the observations of the massive stellar content of R136 suggest 
a physical stellar mass limit near 77T. max * — 150 M Q . 

Furthermore. IWeidner Kroupa ( 20041 ) deduced that the Massey assertion 
would be correct for both cluster masses if the IMF had a slope a.3 ^ 2.8. Such 
a steep slope would make the observed limit consistent with random selection 
from the IMF, and it may be the true power-law index if unresolved multiple 
systems among O stars are corrected for, but this awaits a detailed study. A 
further caveat comes from unresolved multiple systems which would allow an 
?71max*,true Small as ~ W max * /2 if a 3 fa 2.35. 

Similar results were obtained by Don Figer for the Arches cluster. The 
Arches is a star-burst cluster within 30 pc i n projected distanc e from the Galac- 
tic centre. It has a m ass M fa 1 x 10 5 dBosch et alj|200lft . age 2 - 2.5 Myr 
(jNaiarro et al.ll2004l ) and [Fe/H]w ([Najarro et al.ll2004l ). "it is thus a counter- 
part to R136 in that the Arches is metal rich and was born in a very different 

tidal environment to R136. 

Using his HST observations of the Arches, Figer ( 20051 ) performed the same 
analysis as IWeidner &; Kroupal ( 20041 ) did for R136. The Arches appears to 
be dy namically evolved, with sub stantial mass loss through the strong tidal 
forces ( Portegies Zwart et al.ll2002l ) and the stellar mass function with a = 1.9 
is thus flatter than the Salpeter IMF. Using his updated IMF measurement, 
Figer calculated the expected number of stars above 150 M & to be 33, while 
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Figure 4. Number of stars above mass m for R136 with different mass 
estimates (dotted line: Mri- ; 6 = 2.5 x 10 5 M Q , dashed line: Mri 36 = 5 x 
10 4 M Q , (ISelman et al.lll999f). The v ertical solid line marks m = 150 Mq. 
Taken from IWeidner fc Kroupal (|2004h . 



a Salpeter IMF would predict there to be 18 stars. Observing no stars but 
expecting to see 18 has a probability of P = 10 -8 , again strongly suggesting 

"Imax* ~ 150 Mq. 

Given the importance of knowi ng if a finite ph y sical u pper mass limit exists 
and how it varies with metallicity, lOev fc Clarkel (120051 ) studied the massive- 
star content in 9 clusters and OB associations in the MW, the LMC and the 
SMC. They predicted the expected masses of the most massive stars in these 
clusters for different upper mass limits (120, 150, 200, 1000 and 10000 M ). For 
all populations they found that the observed number of massive stars supports 
with high statistical significance the existence of a general upper mass cutoff in 
the range m max * € (120, 200M©)H| 

The general indication thus is that a physical stellar mass limit near 150 Mq 
seems to exist. While biases due to unresolved multiples that may steepen 
the IMF and/or reduce the true maximal mass need to be studied further, the 
absence of variations of m m& ^ with metallicity poses a problem. A constant 
wimax* would only be apparent for a true variation as proposed by the theoretical 
models, if metal-poor environments have a larger stellar multiplicity, the effects 
of which would have to compensate the true increase of m max * with metallicity. 



Maximal stellar mass in clusters Above we have seen that there seems to 
exist a universal physical stellar mass limit. However, an elementary argu- 
ment suggests that star-clusters must also limit the masses of their constituent 



More recent work o n th e physical upper mass limit can be found in iKoenl ((2006); 
iMai'z Apellaniz et all l|2006f l. 
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stars: A pre-star-cluster gas core with a mass Af core can, obviously, not form 
stars with masses m > eM COTe , where e ~ 0.33 is the star-formation efficiency 
(|Lada Ladal I2003T ) . Thus, given a freshly hatched cluster with stellar mass 



M ec i, stars in that cluster cannot surpass masses m max = M ec \, which is the 
identity relation corresponding to a "cluster" consisting of one massive star. As- 
suming the stellar IMF is a continuous density distribution function and that 
clusters are filled with stars distributed according to the stellar IMF, this can 
be generalised by stating that each cluster can have only one most massive star, 

max* 

^(m')dm', (7) 
with 

fmmax 

M ec i(m max ) = / m'£(m')dm' (8) 



-'"Mow 

as a further condition, as above. These two equations need t o be solved numeri- 



cally and give the senii-aiialyticcil relation Tn max — F(M ecX ) (IWeidner fc Kroupal 
2004). It is plotted in Fig. as the thick-solid curve. 



A compilation of clusters from the literature for which the cluster mass and 
the initial mass of the heaviest star can be estimated ( Weidner Sz Kroupa 2006) 



shows that the cluster mass indeed appears to have a limiting influence on the 
stellar mass within it. The observational data are plotted in Fig. \5\ finding 
rather excellent agreement with the semi-analytical description above. 

However, it would be undisputed that a stochastic sampling effect from the 
IMF must be present when stars form. This ca n be mimicked in the com puter 
by performing various Monte-Carlo experiments ([Weidner &: Kroupall200 6j). The 
Monte-Carlo experiments are conducted in three different ways, 

- pure random sampling (random sampling) 

- mass constrained random sampling (constrained sampling) 

- mass constrained random sampling with sorted adding (sorted sampling) 

For the random sampling model, 10 million clusters are randomly taken from 
a cluster distribution with a power-law index of (3m = 2.35 between 12 and 
2.5 x 10 7 stars. The relevant distribution function is the embedded-cluster star- 
number function (ECSNF), 

dN^ocN-P", (9) 

which is the number of clusters containing N G [N',N' + dN') stars. Each 
cluster is then filled with N stars randomly from the standard IMF (eq. [20] 
without a mass limit, or by imposing the physical stellar mass limit, m < 
150 Mq. The stellar masses are added to get the cluster mass, M ec \. For each 
cluster the maximal stellar mass is searched for. For each cluster in a mass bin 
M ec i — AM ec i/2, M ec \ + AM ec j/2 the average m max is calculated, and the set of 
average m max values define the relation 

m max = m^ x (M ecl ). (10) 

For the constrained sampling model, 5 x 10 7 clusters are randomly taken 
from the embedded-cluster mass function (ECMF), 

£ecl(M ecl ) OC Mj, (11) 
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Figure 5. The thick solid line shows the dependence of the mass of the most- 
massive star in a cluster on the cluster mass according to the semi-analytical 
model. The thick dashed line shows the mean maximum stellar mass for sorted 
sampling. The dot-dashed lines are mass-constrained random-sampling results 
with a physical upper mass limit of m max * = 150 M© (thick line) and 10 6 Mq 
(thin line). Pure random sampling models are plotted as dotted lines. The 
thick one is sampled to m max » — 150 Mq while the thin one up to 1O 6 M0. 
The thin solid line shows the identity relation, where a "cluster" consists only 
of one star. The dots with error bars are observed clusters, whi le the triangle 
is a r esult from a s tar-formation simulation w ith an SPH code (Bonnell et al.l 
I2003D . Taken from lWeidner fc Kroupal (|2006f) . 
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between 5M & (the minimal, Taurus- Auriga-type, star-forming "cluster" count- 
ing ~ 15 stars) and 10 7 M & (an approximate maximu m mass for a single s tellar 



population that consists of one metallicity and age, IWeidner et al.ll2004l ) and 
again with f3 = 2.35. Note that (i N « j3 because the ECSNF and the ECMF 
only differ by a nearly-constant average stellar mass. Then stars are taken ran- 
domly from the standard IMF and added until they reach or surpass the respec- 
tive cluster mass, M ec \. Afterwards the clusters are searched for their maximum 
stellar mass. For each cluster in a mass bin M ec \ — AM ec \/2, M ec \ + AM ec ]/2 
the average m max is calculated, and the set of average m max values define the 
relation 

m max = m^Med). (12) 

For the sorted sampling model again 5 x 10 7 clusters are randomly sampled 
from the ECMF (eq.HH) between 5 M© and 10 7 M Q and with = 2.35. However, 
this time the number N of stars which are to populate the cluster is estimated 
from N = M ec i/m av , where m av = 0.36 M & is the average stellar mass for the 
standard IMF (eg. [20]) between 0.01 M & and 150 M & . These stars are added to 

give Me C l iran , 

Median = m, 
N 

such that rrii < m; + i. If M e d 5 ran < M ec \ in this first step, an additional number 
of stars, AN, is picked randomly from the IMF, where AN = {M ec \ — M ran ) / m av 
(we assume m av = constant). Again these stars are added to obtain an improved 
estimate of the desired cluster mass, 



2 



M e d,ran = )^ rUj, < m i+1 . 

N+AN 



This is done such that 2 M e d ira n ~ M ec \ (for details of the method see lWeidner &; Kroupa 
2006]). The procedure is repeated until all clusters from the ECMF are 'filled'. 
They are then also searched for the most massive star in each cluster, as above. 
For each cluster in a mass bin M ec \ — AM ec \/2, M ec \ + AM ec ]/2 the average m max 
is calculated, and the set of average m max values define the relation 

m max = m^(M e d). (13) 

All three relations are plotted in Fig. [5j We noted already that the obser- 
vations follow the semi-analytic relation remarkably well. Furthermore, Fig. [5] 
also suggests that the different Monte-Carlo schemes can be selected for. Thus, 
the sorted-sampling algorithm leads to virtually the same results as the semi- 
analytical relation, and it fits the data very well indeed. The correspondence 
of the sorted-sampling algorithm to the semi-analytical result is not really sur- 
prising, because the algorithm is Monte-Carlo integration of the same problem. 
The constrained-sampling and random-sampling algorithms, on the other hand, 
can be excluded with very high confidence by pe rforming statistical tests o n the 
observational data that a re reported in detail in lWeidner &: Kroupa ( 20061 ) . 



On a historical note, iLarsonl ( 19821 ) had pointed out that more massive and 



dense clouds correlate with the mass of the most massive stars within them 
and he estimated that m max = 0.33 M^^j (masses are in M & ). An updated 
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with the stellar mass 
1.2 M^|^ . Both are flatter than the semi-analytical 



relation was derived by lLarson (|2003l 1 by comparing 
in a few clusters, m ma , 

relat io n, and there f ore do not fit the data in Fig. [5] as well (jWeidner Kroupa 
20061 ). lElmegre"en ( 19831 ) constructed a relation between cluster mass and its 
most massive star based on an assumed equivalence between the luminosity 
of the cluster population and its binding energy, for a Miller- S calo I MF. This 
function is even shallo wer than the one estimated by lLarson (2003) relation. 
Assuming m max * = oo. lElmegreen (2000) solved eqs [7] and [8] above for a single 
Salpeter power-law stell ar IMF finding a rn mm (M er ]) relation quite consistent 
with the data in Fig. [5] (jWeidner Kroupal 120061 ). 



Implications: stellar astrophysics and the formation of star clusters We are 
now in the happier situation that a physical stellar mass limit seems to have 
been found. But the absence of clear variation of this limit with metallicity poses 
a potential problem, although it may be too early to make definite statements. 
Further observational work on many more very young and massive clusters is 
needed to ascertain the findings reported here, and to quantify the multiplicity 
properties of massive stars, as noted above. 

That the sorted-sampling algorithm for making star clusters fits the obser- 
vational maximal-stellar-mass-star-cluster-mass data so well would appear to 
imply that clusters form in an organised fashion. The physical interpretation of 
the algorithm (i.e. of the Monte-Carlo integration) is that as a pre-cluster core 
contracts under self gravity the gas densities increase and local density fluctua- 
tions in the turbulent medium lead to low-mass star formation, perhaps similar 
to what is seen in Taurus-Aurigae. As the contraction proceeds and before 
feedback from young stars begins to disrupt the cloud, star-formation activity 
increases in further density fluctuations with larger amplitudes thereby forming 
more massive stars. The process stops when the most massive star s that have 
just formed supply sufficient feedback energy to disrupt the cloud (jElmegreen 
19831 ). Thus, less -massive pre-cluster cloud-cores would die at a lower maximum 
stellar mass than more massive cores. But in all cases stellar masses are limited, 
m < m max (M ec i) < 

This scenario is nicely co nsistent with t he hy drod ynamic cluster forma- 
tion calculations presented by iBonnell et al. J 2003 ) and [ Bonnell et all (|2004l ) , 
as is reported in mo re detail in IWeidner Kroupal ( 20061 ). We note here that 
Bonnell et al. ( 20041 ) found their theoretical clusters to form hierarchically from 
smaller sub-clusters, and together with continued competitive accretion this 

2 /3 

leads to the relation m max oc M ecl in excellent agreement with the compilation 
of observational data. While this agreement is stunning, the detailed outcome of 
the currently available SPH modelling in terms of stellar multiplicities is not right 
( Goodwin et a.1.1 [2004j : iGoodwin &: Kroupal 120051 ). and feedback that ultimately 
dominates the process of star-formation, given the generally low star-formation 
efficiencies observed in cluster-forming volumes, is not yet incorporated in the 
modelling. This poses a severe computational hurdle because radiation transfer 
requires the correct treatment of atomic processes in remote regions of space 
that are exchanging radiation (eg. iBaes et al.ll2005l 1. 



Caveats Unanswered questions regarding the formation and evolution of mas- 
sive stars remain. There may be stars with m > m max * which implode "in- 
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visibly" after 1 or 2 Myr. The explosion mechanism sensitively depends on the 
prese ntly still rat her uncertain mechanism for shock revival after core collapse 
(e.g. |Jankall200lh . Since such stars would not be apparent in massive clusters 
older than 2 Myr they would not affect the empirical maximal stellar mass, and 
??imax*,true would be unknown at present. 

Furthermore, and as stated already above, stars are often in multiple sys- 
tems. Especially massive stars se em to have a binary fraction of 80% or even 
larger ( Garcia &: Mermilliod 200ll ) and apparently tend to be in binary systems 
with a preferred mass-ratio near unity. Thus, if all O stars would be in equal- 
mass binaries, then m max * t rue ~ "Wx*/ 2 - 

Finally, it is disconcerting that m max * ~ 150 Mq appears to be the same 
for low-metallicity environments ([Fe/H]= —0.5, R136) and metal-rich environ- 
ments ([Fe/H]= 0, Arches), in apparent contradiction to the theoretical values 
(|Stotherdll992T ). Clearly, this issue needs further study. 



4. Low-mass stars 



There are three well-tried approaches to determine f (Afv) in eq. [T] ((Kroupa 
200ld ). The first two are applied to Galactic- field stars, and the third to star 



clusters. The sample of Galactic-field stars close to the Sun is especially im- 
portant because it is the most complete and well-studied stellar sample at our 
disposal. 



4.1. Galactic-field stars 

Galactic-field stars have an average age of about 5 Gyr and represent a mixture 
of many star- formation events. The IMF deduced for these is therefore a time- 
averaged composite IMF. For m ^ 1.3 Mq the composite IMF equals the simple- 
stellar IMF (§E3), and so it is an interesting quantity for at least two reasons: For 
the mass-budget of the Milky- Way disk, and as a bench-mark against which the 
IMFs measured in presently occurring star-formation events can be compared 
to distill possible variations about the mean. 

The first and most straightforward method consists of creating a local 
volume-limited catalogue of stars. Completeness of the modern Jahreiss-Gliese 
Catalogue of Nearby Stars extends to about 25 pc for m ^ 0.6 Mq, trigonometric 
distances having been measured using the Hipparcos satellite, and only to about 
5 pc for less massive stars for which we still rely on ground-based trigonometric 
parallax measurements^. The advantage of the LF, f near (M^), created using 



i Owing to the poor statistical definition of 'I' near (My) for My ^ 10, m^SO.5, it is important 
to increase the sample of nearby stars, but controversy exists as to the maximum distance 
to which the VLMS census is complete. Using spectroscopic parallax it has been suggested 
that the local census of VL MSs is complete to within about 15 % to distances of 8 pc and 
beyond ijReid fc GizislH997l ). However, Malmquist bias allows stars and unre solved binaries to 
enter such a flux- limited sample from much larger distances (|Kroupall2001cD . The increase of 
the number of stars with distance using trigonometric distance measurements shows that the 
nearby sample becomes incomplete for distances larger than 5 pc and for My > 12 (|Jahreissl 
119941 ; [Henry et aDll997l ). The incomple teness in the no rthern stellar census beyond 5 pc and 
withi n 10 pc amounts to about 35 % ( Jao et al.ir2003T ). and recently discovered companions 
(e.g. iDelfosse et all 119991 ; iBeuzit et all 120041 ) to known primaries in the distance range 5 < 
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this catalogue is that virtually all companion stars are known, that it is truly 
distance limited and that the distance measurements are direct. 

The second method is to make deep pencil-beam surveys using photographic 
plates or CCD cameras to extract a few hundred low-mass stars from a hundred- 
thousand stellar and galactic images. This approach leads to larger stellar sam- 
ples, especially so since many lines-of-sight into the Galactic field ranging to 
distances of a few 100 pc to a few kpc are possible. The disadvantage of the LF, 
^phot(My), created using this technique is that the distance measurements are 
indirect relying on photometric parallax. Such surveys are flux limited rather 
than volume limited and pencil-beam surveys which do n ot pass through vir - 
tually the entire stellar disk are prone to Malmquist bias ( Stobie et al. 1989). 



This bias results from a spread of luminosities of stars that have the same colour 
because of their dispersion of metallicities and ages, leading to intrinsically more 
luminous stars entering the flux-limited sample and thus biasing the inferred ab- 
solute luminosities and the inferred stellar space densities. Furthermore, binary 
systems are not resolved in the deep surveys. 

The local, nearby LF and the Malmquist-corrected deep photometric LF are 
displayed in Fig. [6l They differ significantly fo r stars fainter t han My ~ 11.5 
which caused significant controversy in the past (|Kroupalll995al ). That the local 
sample has a spurious but significant over-abundance of VLMSs can be ruled 
out by virtue of the large velocity dispersion in the disk, ~ 30 pc/Myr. Any 
significant overabundance of stars within a sphere with a radius of 30 pc would 
disappear within one Myr, and cannot be created by any physically plausible 
mechanism from a population of stars with stellar ages spanning the age of 
the Galactic disk. The shape of VPphot(-^v) f° r My ft 12 is confirmed by many 
independent photometric surveys. That all of these could be making similar 
mistakes, such as in colour transformations, becomes unlikely on consideration 
of the LFs constructed for completely independent stellar samples, namely star 
clusters (Fig. [7]). 

Eq. [1] shows that any non-linear structure in the MLR is mapped into ob- 
servable structure in the LF, provided the MF does not have compensating 
structure. Such a conspiracy is implausible because the MF is defined through 
the star-formation process, but the MLR is a result of the internal constitution 
of stars. The MLR, its derivative and deviations of models from observational 
data are shown in Figs. [8] and [9J respectively. It is apparent that the slope is 
very small at faint luminosities leading to large uncertainties in the MF near the 
hydrogen burning mass limit. 

The physics underlying the non-linearities of the MLR are due to an in- 
terplay of changing opacities, the internal stellar structure and the equation of 
state of the matter deep inside the stars. Starting at high masses (m ^ fewM ), 
as the mass of a star is reduced H~ opacity becomes increasingly important 
through the short-lived capture of electrons by H-atoms resulting in reduced 
stellar luminosities for intermediate and low-mass stars. The m(My) relation 
becomes less steep in the broad interval 3 < My < 8 leading to the Wielen dip 



d < 12 pc indeed su ggest that t he extende d samp l e may not yet b e complete. Based on the 
work by iReid et al.l |2003f) and iReid et all (|2003T ). iLuhmanl (|2004f ) however argues that the 
incompleteness is only about 15 %. 
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Figure 6. Stellar luminosity functions (LFs) for solar-neighbourhood stars. 
The photometric LF corrected for Malmquist bias and at the mid-plane 
of the Galactic disk ("Jphot) is compared with the nearby LF (\& n ear)- 
The average, g round-based ^phot (dashed histogram, data pre-dating 1995, 
iKroupal Il995al ) is confirmed by Hubble-Space- Telescope (HST) star-count 
data which pass through the entire Galactic d isk and are thus less prone 
to Malmquist bias (solid dots. lZheng et al1l2001f ). The ground-based volume- 
limited trigonometric-parallax sample (dotted histogram) systematically over- 
estimates ^noar due to the Lutz-Kelker bias, thus lying above the im- 
proved estimate provided by the Hip parcos-satellite data (solid histogram, 
iJahreiB fc WielerJll997tlKroupall2001(j) . The depression/plateau near M v = 7 
is the Wielen dip. The maximum near My « 12, Mj s» 9 is the KTG peak. 
The thin dotted histogram at the f aint end indicat es the level of refinement 
provided by recent stellar additions (lKroupall2001cf ) demonstrating that even 
the immediate neighbourhood within 5.2 pc of the Sun probably remains in- 
complete at the faintest stellar luminosities. 
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Figure 7. /-band LFs of stellar systems in fo ur star clusters: globular 
cluster (GC) M15 (|de Marchi fc Parescd Il995al distance m odulus Am = 
m — M = 15.25 mag); GC NGC 6397 (iParesce et al] U995l Am = 12.2); 
young Galactic cluster Pleiades (lHamblv et al.l I1991L Am = 5.48); GC 
47 Tuc (|de Marchi fc Parescd Il995bl Am = 13.35). The dotted histogram 
is ^ P hot(Mj) from the upper panel, transformed to the /-band using the lin- 
ear colour-mag nitude relation M v = 2.9 + 3.4 (V - I) (jKroupa et alJl!993f > 
and * p hot(Mj) = (dMv/dM^^photiMv) (eq.EJ. 
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(Fig. [U]). The m(-My) relation steepens near My = 10 because the formation 
of H2 in the very outermost region of main-sequence stars causes the onset of 
convection up to and above the photo-sphere leading to a flattening of the tem- 
perature gradient and therefore to a larger effective temperature as opposed to 
to an artificial case without convection but the same central temperature. This 
leads to brighter luminosities and full convection for m < 0.35 Mq. The modern 
Delfosse data beautifully confirm the steepening in the interval 10 < My < 13 
predicted in 1990, and the dotted MLR demonstrates the effects of suppress- 
ing the formation of the H2 molecule by lowering it's dissociation energy from 



4.48 eV to 1 eV dKroupa et alJll990h . The m{M v ) relation flattens again for 



My > 14, m < 0.2 M Q as degeneracy in the stellar core becomes increasingly im- 
portant for smaller masse s limiting further contraction ( Havashi Nakandl 19631 : 
Chabrier Sz Baraffd 1997). Therefore, owing to the changing conditions within 
the stars with changing mass, a pronounced local maximum in —dm/dMy re- 
sults at My ~ 11.5, postulated in 1990 to be the origin of the maximum in ^phot 



near My = 12 dKroupa et alJll990l ) 



The implication that the LFs of all stellar populations should show such a 
feature, although re alistic metallicity-d ependent stellar models were not avail- 
able yet, was noted ( Kroupa et al. 19931 ). The subsequent finding that all known 
stellar populations have such a maximum in the LF (Figs. [6] and [7J constitutes 
one of the most impressive achievements of stellar- structure theory. Different 
theoretical m(My) relations have the maximum in —dm/dMy at different My, 
suggesting the possibility of testing stellar structure th eory near the critical 
mass m gj 0.35 M^, where stars become fully convective (jKroupa Toutlll997l ; 
Brocato et al.l fl998l ) . But since the IMF also defines the LF the shape and loca- 



tion cannot be unambiguously used for this purpose unless it is postulated that 
the IMF is invariant. 

A st udy of the position of th e ma ximum in the /-b a nd LF has been under- 
taken by Ivon Hippel et al.1 (ll996h and iKroupa fc Tout! (| 19971 ) finding that the 
observed position of the maximum shifts to brighter magnitude with decreasing 
metallicity, as expected from theory (Figs. [TU1 and ITT|) . 

In addition to the non-linearities in the m(Mp) relation, unresolved mul- 
tiple systems affect the MF derived from the photometric LF, in particular 
since no stellar population is known to exist that has a binary proportion 
smaller than 50 per cent, apart possibly from dynamically highly evolved glob- 
ular clusters. Suppose an observer sees 100 systems. Of these 40, 15 and 
5 are binary, triple and quadruple, respectively, these being realistic propor- 
tions. There are thus 85 companion stars which the observer is not aware 
of if none of the multiple systems are resolved. Since the distribution of sec- 
ondary masses for a gi yen primary mass is not un iform but typically increases 
with decreasing mass ( Malkov Zinneckerl l200ll ). the bias is such that low- 
mass stars are si gnificantly underrepresented in any survey that does not de- 
tect companions (Kroupa et al.lll99ll ; lHoltzman et alJll99St iLuhman et ai1ll998l : 
Malkov fc Zinneckerll200l[ r Also, if the companion(s) are bright enough to affect 



the system luminosity noticeably then the estimated photometric distance will 
be too small, artificially enhancing inferred space densities which are, however, 
mostly compensated for by the larger distances sampled by binary systems in a 
flux-limited survey, together with an exponential density fall-off perpendicular 
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Figure 8. The mass-luminosity relation (MLR, upper panel) and its deriva- 
tive (lower panel) for late-type stars. Uppe r panel: The most r ecent observa- 
tional data (sol id triangles and open circles. iDelfosse et al.ll2000t open squ ares, 
lAndersenlll991h are compared wit h the empirical ML R of IScalo (1986) and 
the semi- empirical KTG93-MLR (jKroupa et alJ l"l993). The under- luminous 
data points 1,2 are GJ2069 Aa,b and 3,4 are G1791.2A,B. All are proba- 
bly metal-rich by ~0-5 d ex (|Delfosse et al.ll2000D. Recent theoretical MLRs 
from lBaraffe et al.l (|1998l ) (BCA H98) andlSiess et all (|2000D (SDFOO) are also 
shown. The observational data (|Andersenlll991h show that logio [to (My)] is 
approximately linear for m > 2 M Q . Lower panel: The derivatives of the same 
relations plotted in the upper panel are compared with ^phot from Fig.0 Note 
the good agreement between the location, amplitude and width of the KTG 
peak in the LF and the extremum in dm/dMy- 
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Figure 9. Deviations of the MLRs (dm = 77i b s — m(Mv)) from the em- 
pirical data with errors e m shown in Fig. [8] in Mq (upper panel) and in 
percent (lower panel with uncertainties — (m(Mv) / m^ ha ) X e m ). Reduced xl 
(v = 26 for 31 data points, ignoring the four outliers) values indicate the 
formal goodness-of-fit. Formally, none of the MLRs available is an accept- 
able model for the data. This is not alarming though, because the models 
are for a single-metallicity, single-age population while the data span a range 
of metallicities and ages typical for the solar-neighbourhood stellar popula- 
tion, as signified b y dm ^> e m in most cases. The \t values confirm that the 
BCAH98 models fearaffe et al.l [l998h and the semi-empirical KTG93 MLR 
()Kroupa et al.lll993h provide the best-matching MLRs. Note that the KTG93 
MLR was derived from mass- luminosity data prior to 1980, but by using the 
shape of the peak in v I / phot(-^y) as an additional constraint the constructed 
MLR became robust. The lower panel demonstrates that the deviations of 
observational data from the model MLRs are typically much smaller than 
30 per cent, excluding the putatively metal-rich stars (1-4). 
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Figure 10. The position of the maximum in —d m/dMi as a func t ion o f 
metallicity of the theoretical mass-luminosity data of iKroupa &: Toutl (|1997[ ) 
is shown as solid d ots. The open squares r epres ent bounds by the stellar- 
structure models of iD'Antona fc Mazzite lli (1996), and the open circles are 
observational constraints for different populations (e.g. SN for the compos- 
ite solar-neighbourhood population, Pleiades for the simple population of an 
intermediate-age cluster). Thick circles are more certain than the thin circles, 
and for the Pleiades a sequence of positions of the LF-maximum is given, 
from top to bottom, with the following combinations of (distance modulus, 
age): (5.5, 70 Myr), (5.5, 120 Myr), (5.5, main se quence), (6, 70 Myr) (6, 
120 Myr), (6, main sequence). For more details see IKroupa fc Toutl (fl99l . 
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Figure 11. Similar to Fig. Qj}]but from lvon Hippel et al.l (|1996l ). their fig. 5: 
The absolute /-band-equivalent magnitude of the maximum in the LF as a 
function of metallicity for different populations. The solid and dashed lines 
are loci of constant mass ( 0.2, 0.231, 0.3 Mm ) accor ding to theoretical stellar 
structure calculations. See Ivon Hinnel et all (fl99rl for more details. 



to the Galactic disk ( Kroupa 2001c! ) . A faint companion will also be missed if 
the system is formally resolved but the companion lies below the flux limit of 
the survey. 

Comprehensive star-count analysis of the solar neighbourhood that incor- 
porate unresolved binary systems, metallicity and age spreads and the density 
fall-off perpendicular to the Galactic disk show that the IMF can be approx- 
imated by a two-part power-law with a% = 1.3 ± 0.7,0.08 < tu/Mq < 0.5, 
a? = 2.2, 0.5 < tu/Mq < 1, a result obtained for two different MLRs ([Kroupa 



200ld h The index 

«2 is constrained tightly owing to the well-constrained ^ r 



the well-constrained empirical MLR in this mass range and because unresolved 
binary systems do not significantly affect the solar-neighbourhood LF in this 
mass range because primaries with m ^ 1 M & are rare and are not sampled. 

Fig. [12] demonstrates models of the single-star and system LFs for the 
KTG93 MLR shown in Fig. [5J The significant difference between the single- 
star and system LFs is evident, being most of the explanation of the disputed^ 



9 The discrepancy evident in Fig. 1121 between the nearby LF, f near , and the photometric LF, 
typhot, invoked a significant dispute as to the nature of this discrepancy. On the one hand 
(Krouoa 1995a) the difference is thought to be due to unseen companions in the deep but low- 
resolution surveys used to construct *& p hot, with the possibility that photometric calibration 
for VLMSs may remain proble matical so that the exact shape of Vfphot for My <^ 14 is probably 
uncertain. On the other hand (Reid & Gizis 1997) the difference is thought to come from non- 
linearities in the V — I, My colour-magnitude relation used for photometric parallax. Taking 
into account such structure it can be shown that the photometric surveys underestimate stellar 
space densities so that 'I'phot moves closer to the extended estimate of 'I'ncar using a sample of 
stars within 8 pc or further. While this is an important point, the extended ^near is incomplete 
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Figure 12. Comparison of the model field luminosity function (curves) of 
a single-metallicity and single-age population that is without measurement 
errors, with observations in the photometric V-band (a co mparison of th e 
corresponding LFs in bolometric magnitudes can be found in lKroupalll995el ). 
The model assumes the standard stellar IMF, eq. [20] below. The model sin- 
gle star luminosity function is normalised to the nearby luminosity function 
at My ~ 10, Mboi ~ 9, giving the normalisation constant in the MF k, and 
the plot shows fc'I'mod.sing (long dashed curve), k vE r m od,sys(* = 0) (dotted 
curve, without pre-main sequence brightening) and k ^ mo d,sys(t = 1 Gyr) 
(solid curve) . Note that the solid curve is the luminosity function for a realis- 
tic model of the Galactic field population of systems consisting of 48 per cent 
binaries which have a period distribution consistent with the empirical G-, 
K-, and M-dwarf peri od distribution, the mass ra tio distributions for G-dwarf 
systems as observed (iDuauennov k, Mayor! 199 j|) . and the overall mass-ratio 
distribution given by fig. 2 in lKroupa et al.1 ( 2003T ). where a concise description 
of the "standard star-formation model" can be found. The observed nearby 
stellar luminosity func tion, 1^^, wh ich is not corrected for Malmquist-type 
bias (tables 2 and 8 in lKroupal (|1995af ) and which is smoothed by using larger 
bin widths at the faint end, as detailed in section 4 of that paper, is plot- 
ted as the solid-line histogram. The filled circles represent the best-estimate 
Malmquist corrected photom etric luminos i ty fun ction, ^phot, Fig. [6] By cor- 
recting for Malmquist bias (|Stobie et al.l 119891 ) the LF becomes th at of a 
single-age, single-metallicity population. Taken from iKroupal (|1995eh . 
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discrepancy between the observed 'I' near and ^phot- Note though that the ob- 
served photometric LF contains triple and quadruple systems that are not ac- 
counted for by the model. Note also that the photometric LF has been corrected 
for Malmquist bias and so constitutes the system LF in which the broadening 
due to a metallicity and age spread and photometric errors has been largely 
removed. It is therefore directly comparable to the model system LF, and both 
indeed show a very similar KTG peak. The observed nearby LF, on the other 
hand, has not been corrected for the metallicity and age spread nor for trigono- 
metric distance errors, and so it appears broadened. The model single-star LF, 
in contrast, does not, by construction, incorporate these, and thus appears with 
a more pronounced maximum. Such observation al effects can be inc orporated 
rather easily into full-scale star-count modelling ( Kroupa et al.|[l993T ). The de- 
viation of the model system LF from the observed photometric LF for My <^ 14 
may indicate a change of the pairing properties of the VLMS or BD population 

Since the nearby LF is badly defined statistically for My ~ 13, the resulting 
model, such as shown in Fig. [T2l is a prediction of the true, single-star LF 
that should become apparent once the immediate solar-neighbourhood sample 
has been enl arged significantly t hrough the planned space-based astrometric 
survey Gaia ( Gilmore et al. 19981 ). followed by an intensive follow-up imaging 



and radial-velocity obs erving program me scrutinising every nearby candidate 
for unseen companions ( Kroupall200ld ). Despite such a monumental effort, the 
structure in will be smeared out due to the metallicity and age spread of 

the local stellar sample, a factor to to be considered seriously. 



4.2. Star clusters 

Star clusters offer populations that are co-eval, equi-distant and that have the 
same chemical composition, but as a compensation of these advantages the ex- 
traction of faint cluster members is very arduous because of contamination by 
the background Galactic-field population. The first step is to obtain photometry 
of everything stellar in the vicinity of a cluster and to select only those stars that 
lie near one or a range of isochrones, taking into account that unresolved bina- 
ries are brighter than single stars. The next step is to measure proper motions 
and radial velocities of all candidates to select only those high-probability mem- 
bers that have coinciding space motion with a dispersion consistent with the a 
priori unknown but estimated internal kinematics of the cluster. Since nearby 
clusters for which proper-motion measurements are possible appear large on the 
sky, the observational effort is horrendous. For clusters such as globulars that 
are isolated the second step can be omitted, but in dense clusters stars missed 
due to crowding need to be corrected for. The stellar LFs in clusters turn out 
to have the same general shape as ^phot (Fig. [7|), with the maximum being 



(see footnote [8]) and theoretical colour-magnitude relations do not have the required degree 
of non-linearity. The observational c olour-magnitude d ata also do not conclusively suggest 
a feature with the required strength ijBaraffe et al.lfl998f ). Furthermore, ^ p hot agrees almost 
perfectly with the LFs measured for star clusters of solar and population I metallicity (Fig. [TJl so 
that it appears unlikely that non-linearities in the colour-magnitude relation are the dominant 
source of the discrepancy. 
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slightly off-set depending on the metallicity of the population (Figs. [TU] and [HQ. 
A 100 Myr isochrone (the age of the Pleiades) is also plotted in Fig. [8] to empha- 
sise that for young clusters additional structure (in this case another maximum 
near M y = 8 in the LF is e xpected via eq. Q] . This is verified for the Pleiades 
cluster ( Belikov et al. 19981 ). a nd is due to stars with m < 0.6 Mq not having 



reached the main-sequence yet (jChabrier k, Bara ffc 2000) 



LFs for star clusters are, like ^photi system LFs because binary systems are 
not resolved in the typical star-count survey. The binary-star population evolves 
due to encounters, and after a few initial crossing times only those binary systems 
survive that have a binding energy larger than the typical kinetic energy of stars 
in the cluster. A further complication with cluster LFs is that star clusters 
preferentially loose single low-mass stars across the tidal boundary as a result 
of ever-continuing re-distribution of energy during encounters while the retained 
population has an increasing binary proportion and increasing average stellar 
mass. The global PDMF thus flattens with time with a rate proportional to 
the fraction of the cluster lifetime and, for highly evolved initially rich open 
clusters, it evolves towards a delta function near the turnoff mass, the mass-loss 
rate being a function of Galactocentric distance. This is a major problem for 
aged open clust ers (initially N < 10 4 stars) with life-times of only a few 100 Myr 
(jKrouDdll995ch . 



These processes are now well quantified ( Kroupall2001al ; |Portegies Zwart et~aT 



l200ll ) , and Fig. [13] shows that a dynamically very evolved cluster such as the 
Hyades has been depleted significantly in low-mass stars. Even so, the binary- 
star correction that needs to be applied to the LF in order to arrive at the 
single-star present-day LF is sig nificantB. 

A computationally challenging investigation of the systemat ic changes of the 

MF in evolving clusters of different masses has been published bv lBaumgardt &: Makinol 
(|2003l ). Baumgardt & Makino quantify the depletion of the clusters of low-mass 
stars through stellar-dynamical processes with a hitherto not available quality 
and conclusively show that highly evolved clusters have a very substantially 
skewed present-day MF (Fig. I14|) . If the cluster ages are expressed in fractions, 
Tf, of the overall cluster lifetime, which depends on the initial cluster mass, its 
concentration and orbit, then different clusters on different orbits lead to vir- 
tually the same present-day MFs at the same Tf. Their results were obtained 
for clusters that are initially in dynamical equilibrium and that do not con- 
tain binary stars (these are computationally highly demanding), so that future 
analysis, including initially nqn-vir ialised clusters and a high primordial binary 
fraction ( Goodwin &: Kroupall2005l ). will be required to further refine these ex- 



cellent results. 

The first realistic calculations of the formation of an open star cluster such 
as the Pleiades demonstrate that the binary properties of stars remaining in 
the cluster are comparable to those observed even if all stars initially f orm i n 
binary systems with Taurus- Auriga- like properties (jKroupa et al.l 12001 . 20031 ) . 



They also demonstrate the complex and counter-intuitive interplay between the 
initial concentration, mass segregation at the time of residual gas expulsion, and 



10 Note that a system can be a binary star or a single star, and that the single-star LF is the LF 
obtained by counting all stars independently of whether they are in systems or single. 
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Figure 13. Models of the if-band single-star and system LFs in an en- 
semble of 20 dynamically highly evolved clusters (thin and thick histograms, 
respectively). Each cluster model consists initially of 200 binaries with a half- 
mass radius of 0.85 pc, and the LFs are shown at an age of 480 Myr (44 initial 
crossing times) and count all stars and systems within the central sphere with 
a radius of 2 pc. The clusters are random renditions from the same parent 
distributions (binary-star orbital parameters, IMF, stellar positions and ve- 
locities) and are initially in dynamical equilibrium. The upper dashed curve 
is the initial single-star LF (KTG93 MLR and standard IMF, eg . 1201 below) 
and the solid curve is the model Galactic-field LF of systems. This is an accu- 
rate representation of the Galactic-field population in terms of the IMF and 
mixture of single and binary stars, and is derived by stars forming in clusters 
such as shown here that dissolve with time. Both of these LFs are identical to 
the ones shown in Fig. [T2l The dotted curve is the initial system LF (100 % 
binaries). From iKroupal (|1995d) . 
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Figure 14. Present-day mass functions in a King-model cluster with concen- 
tration Wo = 5 on a circular orbit about the MW centre with radius 8.5 kpc. 
Shown are the MFs of all bound stars at ages corresponding to rt =0 %, 
30 %, 60 % and 90 % of the cluster life-time (from top to bottom). For 
each age the solid line represents one computation with 1.28 x 10 5 stars, the 
dashed lines show the sum of four clusters each with 8000 stars (scaled to the 
same number of stars as the massive computation). Results for other circular 
and e ccentric orbits and cluster con centrations are virtually indistinguishable. 
From lBaumgardt fc Makinol (|2003t) . 
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Figure 15. The observed MF in the Ple iades cluster. Left p anel: The sym- 
bols are observational data (for details see lMoraux et al1l2004T ) and the curve 
is a log-normal fit. Right panel: The solid curve is the same log- normal fit. 
Theor etical system MFs fo r two initial models of the Pleiades cluster accord- 
ing to iKroupa et al.l (|200l[ ) are plotted at an age of 100 Myr. These models 
assume the young cluster to be deeply embedded in gas, with a star-formation 
efficiency of 33 %, a gas-expulsion time-scale shorter than the crossing time 
and to contain 10 4 stars and BDs. Model A (dashed histogram) has an ini- 
tial central number density pc — 10 4 8 stars/pc 3 , while model B (dotted 
histogram) has pc — 10 5 8 stars/pc 3 . The embedded phase lasts 0.6 Myr, 
and during this time mass segregation develops in the initially denser cluster 
model B. Note that these models are not a fit but a p rediction of the Pleia des 
MF, assuming it had a standard IMF (eq.l2"0]). From lMoraux et all (|2004f l. 



the final ratio of the number of BDs to stars (Fig. [T5l) . Thus, the modelling by 
Kroupa et all (|200lh shows that an initially denser cluster evolves to significant 



mass segregation when the gas explosively leaves the system. Contrary to naive 
expectation, according to which a mass-segregated cluster should loose more 
of its least massive members during expansion after gas expulsion, the ensuing 
violent relaxation of the cluster retains more free-floating BDs than the less- 
dense model. This comes about because BDs are split from the stellar binaries 
more efficiently in the denser cluster. These issues remain an active area of 
research, because at least two changes need to be made to t he modelling: on the 
one hand, BDs need to b e treated as an extra population ( Kroupa et al.l 2003; 



Kroupa &: Bouvierl 120031 ) so that the above discovered free-floating BDs that 



result from the disruption of star-BD binaries will not be available, and on the 
other hand some observations suggest that star clusters may form highly mass- 
segregated (but see Fig. [T7|). The mass-dependent loss of stars thus definitely 
remains an issue to be studied. 

From the above work it is now established that even clusters as young as the 
Pleiades are significantly evolved becau se clusters of a ll masses form from highly 
concentrated embedded morphologies ( Kroupall2005l ). Also, the low-mass stars 
in clusters as young as the Pleiades or M35 (Fig. 1161 below) have not yet reached 
the main sequence, so that pre-main sequence stellar-evolution calculations have 
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to be resorted to when transforming measured luminosities to stellar masses via 
the MLR. 

For ages younger than a few Myr this becomes a serious problem: Classi- 
cal pre-main sequence theory, which assumes hydrostatic contraction of spher- 
ical non-, sometimes slowly-rotating stars from idealised initial states breaks 
down because of the overlap with the star formation processes that defies de- 
tailed treatment. Stars this young remember their accretion history, invali- 
dating the application of classical pre-main sequence stellar evolution tracks, a 
point made explicitly clear by the excellent work of Wuch terl and collaborators 
(jWuchterl Klessenll200ll : IWuchterl fc Tscharnuted 12003). and are m any case 
rotating rapidly and non-spherical. Such realistic pre-main sequence tracks are 
not available yet. 

Nevertheless, owing to the lack of an alternative, research of the IMF in very 
young clusters has to resort to spectroscopic classification of individual stars to 
place them on a theoretical isochrone of existing classical theory to estimate 
masses (e.g. Mever et al.ll2000l ; lLuhmanll2004l ; iBarrado v Navascues et al1l2004l ; 
ISlesnick et al.1 l2004f ). In such cases the age spread becomes comparable to the 
age of the cluster (eq. [3]). Binary systems are mostly not resolved. 

A few results are shown in Figs . 1161 and 1171 While th e usual argument is 
for an invariant IMF (|Kroupal l2001al . 120021: IChabrier|[2003al ). as is apparent for 
most population II stars (e.g. fig. 5 in ChabrieiT 2003a ) . Fig. 1161 shows that some 
appreciable differences in measured MFs are evident. The M35 MF appears to 
be highly deficient in low-mass stars. This clearly needs further study because 
M35 and the Pleiades appear to be otherwise fairly similar in terms of age, 
metallicity (M35 is somewhat less metal-rich than the Pleiades) and the size of 
the survey volume. 

Taking the ONC as the best-studied example of a very young and nearby 
rich cluster (age 1 Myr, distance sa 450 pc; N « 5000 — 10000 stars and 
BDs: iHillenbrand Carpenter! l2000l : iLuhman et aD l2000l : iMuench et aD l2000l : 
Kroupa 20001 : Slesnick et al.ll2004l ). Fig. [TBI shows how the shape of the deduced 



IMF varies with improving (but still classical) pre-main sequence contraction 
tracks. This demonstrates that any sub-structure cannot, at present, be relied 
upon to reflect possible underlying physical mechanisms of star formation. 

For the much more massive and long-lived globular clusters (N ^ 10 5 stars) 
theoretical stellar-dynamical work shows that the MF measured for stars near 
the cluster's half-mass radius is approximately similar to the global PDMF, while 
inwards and outwards of this radius the MF is flatter (small er ai ) and steeper 
(larg er a%), respectively, owing to dynamical mass segregation (jVesperini fc Heggie 
19970 . However, mass loss from the cluster flattens the global PDMF such that 
it no longer resembles the IMF anywhere (Fi g. [14J , for which evidence has been 



found in some cases ( Piotto Zoccali 1999I ). The MFs measured for globular 
clusters must therefore generally be flatter than the IMF, which is indeed born- 
out by observations (Fig. [20] below). However, again the story is by no means 
straightforward, because globular clusters have signific antly smaller binary frac- 
tions than population II clusters ( Ivanova et al.ll2005l ). and so the binary-star 
correct ion is smaller fo r globular cluster MFs. Therefore, and as already pointed 
out bv lKroupal (|2001bh . it appears quiet realistically possible that population II 
IMFs were in fact flatter (smaller a\) than population I IMFs, as would be 
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qualitatively expected from simple fragmentation theory (§ [T7I below). Clearly, 
this issue needs detailed investigation which, however, is computationally ex- 
tremely highly demanding, requiring the use of state-of-the art iV-body codes 
and special-purpose hardware. 

The IMFs obtained for nearby and young star clusters as well as for globular 
clusters are summarised in Table [TJ 



5. Very low-mass stars and brown dwarfs 

These are stars near to the hydrogen-burning mass limit (VLMS) or below (BD), 
in which case they are not massive enough to achieve sufficiently high central 
pressures and temperatures to stabilise against continued contraction by burning 
H and thus indefinitely cool to unobservable luminosities and temperatures. Ob- 
servationally it is very difficult to distinguish between VLMSs and BDs, because 
a sufficiently young BD may have colours and spectral featur es corresponding to 
a VLM S. BDs were studied as theoretical objects in 1963 bv lHavashi &: Nakand 
(1963), who performed the first truly self-consistent estimate of the minimum 
hydrogen burning mass limit by computing the luminosity at the surface and 
the energy release rate by nuclear burning. Modern theory of the evolution and 
internal constitution of BDs has advanced considerably owing to the inclusion 
of an improved equation of state and realistic model-atmospheres that take into 
account absorption by many molecular species as well as dust allowing the identi- 



fication of characteristic photometric signatures (jChabrier k, Baraf fc 2000). The 
first BDs were detected in 1995, and since the n they have been found in the so- 
lar neighbourhood and in young star clusters (Basri 2000) allowing increasingly 



sophisticated estimates of their mass distribution. 

For the solar neighbourhood, near-infrared large-scale surveys have now 
identified many dozens of BDs probably closer than 25 pc ( Allen et al. 20051 ) . 



Since these objects do not have reliable distance measurements an ambiguity 
exists between their ages and distances, and only statistical analysis that relies 
on an assumed star -formation hist ory for the solar neighbourhood can presently 
constrain the IMF ( Chabrierll2002l ) . finding a 60 % confidence interva l qq = 0.3± 



0.6 fo r 0.04- 0.08 M approximately for the Galactic-field BD IMF ([Allen et al.1 
20051 1. 



Surveys of young star clusters have also discovered BDs by finding objects 
that extend the colour-magnitude relation towards the faint locus while being 
kinematical members. Given the great difficulty of this endeavour, only a few 
clusters now possess constraints on the MF. The Pleiades star cluster has proven 
especially useful, given its proximity (d w 127 pc) and young age (t c \ 100 Myr). 
Results indicate ao ~ 0.5 — 0.6 (Table [I]). Estimates for other clu sters (ONC, 
( j Ori , IC 348; Table [1]) also indicate ao~0.8. In their table 1, I Allen et al.l 
(2005) summarise the available measurements for 11 populations finding that 



ao ~ — 1 
5.1. BD binaries 

The above estimates of the IMF suffer under the same bias affecting stars 



name ly unseen co mpanions. BD-BD binary systems are known to exist ( Basri 



2000|), in the field (|Bouv et al.ll2003l : IClose et al.ll2003h and in clusters (iMartin et al 
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Figure 16. Left panel: T he measured system mass fu nctions (lm = 
log 10 (m/M©)) in the ONC ([Hillenbrand &; Carpenter] l2000h : optical data, 
r < 2.5 pc, th < 2 M yr, [Fe/H]= -0.02 (lEsteban eTai]ll99l in the Pleiades 
(Ha mblv et allfl999l): r < 6 7 pc r c i « 100 Myr, [Fe/H]= +0.01, and in M35 
(jBarrado v Navascues et~alll2001l ): r < 4.1 pc, r d w 160 Myr, [Fe/H]= -0.21, 
where r is the approximate radius of the survey around the cluster centre and 
t c \ the nuclear age. The strong decrease of the M35 MF below m « 0.5 M 
remains present despite using different MLRs (e.g. DM97, as in the right 
panel). None of these MFs are corrected for unresolved binary systems. The 
standard single-star IMF (eq. |2"0")) is plotted as the three straight lines. Right 
panel: The shape of the ONC IMF differs significantly for m < 0.22 M© if 
different pre-main sequence evolution tracks, and thus essentially different the- 
oretical MLRs, a re employed (DM stands for tr acks calculated by D'Antona 
& Mazzitelli, see iHillenbrand fc Carpente"rl [20"00' for details.) 
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Figure 17. Left panel: Mass segregation is very pronounced in the ONC, as 
is evident by comparing the M F for all stars within r = 2.5 p c with the MF 
for all stars with r < 0.35 pc ([Hillenbrand fc Carpente^l2000L HCOO) (lm = 
log-in (m/M®)). For both samples the reddening Ay < 2.5 mag ((Hillenbrand 
119971 H97, is for an optical and spectroscopic survey, whereas HCOO is a near- 
infrared survey). Right panel: The ratio, < 2.5 pc)/^^ < 0.35 pc), 
of the MFs shown in the left panel increases significantly with decreasing 
mass, demonstrating the significant depletion of low-mass sta rs in the cen- 
tral r egion of the ONC. Stellar-dynamical models of the ONC (|Kroupa et al.l 
120011 ) approximately reproduce this trend at an age of 2 Myr for the stan- 
dard IMF (eq. [201 whereby the system masses of surviving binary systems 
are counted instead of the individual stars many of which are in unre- 
solved binaries) even if no initial mass segregation is assumed (at t = 0, 
< 2.5pc)/£ l(V < 0-35 pc) = con stant). The model snapshots shown are 
from model B in lKroupa et alj (|2001h under the assumption that prior to gas- 
expulsion, 2-1.4 Myr ago, the central stellar density was pc = 10 5 8 stars/pc 3 . 
The dotted lines are eye-ball fits to the plotted data. 
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2003). Their frequency is not yet very well constrained since detailed scrutiny of 



individual objects is time-intensive on large telescopes. But the results show con- 
clusively that the semi-major axis distribution of VLMSs an d BDs is mu c h mor e 
compact than tha t of M dwarfs , K dwarfs and G dwar f s. iBouv et al.l (2003); 
IClose et al.1 (|2003l ) : iMartm et al.l (|2003l ): lPhan-Bao et al.l (|2005|) all find that BD 
binaries with semi-major axis a ^ 15 AU are very rare. Using Monte-Carlo ex- 
periments on publ i shed m ultiple-epoch radial-velocity data of VLMSs and BDs, 
Maxted &: Jeffries] ( 20051 ) deduce an overall binary fraction between 32 and 45 % 



with a semi-major-axis distribution that peaks near 4 AU and is truncated at 
about 20 AU. In the Pleiades cluster where their offset in t he colour-magnitude 
diagram from the single-BD locus makes them conspicuous, iPinfield et ail ([20031 ) 
find the BD binary fraction may be as high as 60 %. 

The truncated semi-major- axis distribution of BDs may be a result of 



binary -disruption in dense clusters of an initial stellar-like distribution. [Kroupa et al 



(2003) test this notion by setting-up the hypothesis that BDs form according to 
the same pairing rules as stars, the star-like hypothesis. This ought to be true 
since objects with masses 0.04 — 0.072 M Q should not have very different pair- 
ing rules than stars that span a much larger range of masses (0.1 — IMq) but 
show vir tually the same period-distribution function (the M-, K- and G-dwarf 
samples, iFischer & Marcvlll992l ; iMavor et ailll992l ; iDuquennov k Mavor|[l99ll . 



respectively). Thus, the hypothesis is motivated by observed orbital distribution 
functions of stellar binaries not being sensitive to the primary mass, which must 
come about if the overall physics of the formation problem is similar. Further 
arguments for a star-like origin of BDs comes from the detection of accretion 
onto and disks about very young BDs, and that the BDs and s tars in Taurus- 
Aurig a have indistinguishable spatial and velocity distributions ([White Basri 



200. 



iKroupa et al.l (|200ll . [20031 ) perform iV-body calculations of ONC- and Taurus 
Auriga-like stellar aggregates to predict the semi-major-axis distribution func- 
tions of BD-BD, star-BD and star-star binaries. These calculations demon- 
strate that the binary proportion among BDs is smaller than among low-mass 
stars after a few crossing times, owing to their weaker binding energies. The 
distribution of separations, however, extends to similar distances as for stel- 
lar systems (a rj 10 3 AU), disagreeing completely with the observed distribu- 
tion. The star-like hypothesis thus predicts far too many wide BD-BD bina- 
ries. This can also be seen from the distribution of binding energies of real 
BD binaries. It is very different to that of stars by having a low-energy cut- 
off, BD £ b in,cut ~ -lO-°- 9 M (pc/Myr) 2 , that is much higher than that of the 
M dwarfs, M -E b in,cut ~ -10~ 3 M (pc/Myr) 2 (Fig. [TS]). This is a very strong in- 
dicator for some fundamental difference in the dynamical history of BDs which 
sets a different energy scale than for stars. 

Furthermore, the theoretical distributions contain a substantial number of 
star-BD pairs, which also disagrees with the obser vational constrains that have 



found very few BD companions to nearby stars (jBasril 120001 ; iPhan-Bao et al 



2005). Basically, if BDs form exactly like stars, then the number of star-BD bi- 



naries would be significantly larger than is observed, since for example G-dwarfs 
prefer to pair with M-dwarfs. The observed general a bsence of BD companions 
is referred to as the BD desert ( Zucker Sz Mazeh 200 ll ). since stellar companions 
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Figure 18. The distribution of binding energies, E^in = -Gmi m2/(2a), 
of BDs (solid line) compared to those of M dwarfs (dot-dashed line). The 
BD d istribution i s comp uted by adopting the semi-major-axis, a, distribution 
from lClose et all (|2003l ) and choosing 10 7 BD masses, m, S (0.04 — 0.1 M Q ), 
from a power-law MF with ctQ — 0.3. The M d warf energy distr i bution 
is computed by assuming the a-distribution from Fisc her &: Marcvl (|1992f ) 
(which is practically identical to that of G dwarfs) and choosing 10 7 masses, 
rrii £ (0.1 — 0.5 M Q ), from a power-law MF with a% = 1.3. The BD distribu- 
tion is incomplete to the right of the soli d line because tight BD pairs could 
not be discovered by the available surveys (|Bouv et alj|2003t IClose et al]|2003t 
iMaxted fc Jeffrieil2005h . 




and planets are found at such separations ( Halbwachs et alj|200d : IVogt et al 
2002). A few very wide systems with B P companions can form dur ing the final 
stages of dissolution of a small cluster (jde La Fuente Marcos ! 19981). and thre e 
such common proper-motion pairs have perhaps been found (|Gizis et al.ll200ll ). 

Finally, the star-like hypothesis predicts far too few star-st ar binaries in 
Tauru s-Auriga, where binary disruption has not been active. iKroupa et al.l 
(2003) thus conclude that the observed BD population is incompatible with 
their hypothesis, such that BDs need to be treated as a separate, or extra, 
population. 

5.2. The number of BDs per star 

Briceho et al. ( 20021 ) report that Taurus- Auriga appears to form significantly 



fewer BDs per star than the ONC. Both systems are very different physically 
but have similar ages of about 1 Myr. This finding was interpreted to be the first 
possible direct evidence of a variable IMF, being nicely consistent qualitatively 
with the Jean-mass, 



Mj oc p- 1/2 T 3/2 , 



(14) 



being larger in Taurus-Auriga than in the ONC because its gas density, p, is 
smaller by one-two orders of magnitude, while the temperatures, T, are similar 

to within a factor of few. 

Given this exciting finding. [Kroupa et al.l ( 20031 ) computed iV-body models 
of the stellar aggregates in Taurus- Auriga in order to investigate the hypothesis 
that BDs form star-like. They find that the same initial number of BDs per 
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star in Taurus- Auriga and in the ONC actually leads to different observed ra- 
tios because BD-BD and star-BD binaries are disrupted more efficiently in the 
ONC; the observer thus sees many more BDs there than in the comparatively 
dynamically unevolved Ta urus- Auriga gr o ups. But, as already noted above, this 
approach fails because the Kroupa et al. ( 20031 ) model has too many wide BD 
BD binaries, and al s o it pr edicts too many star-BD binaries. Given this problem, 
iKroupa Sz Bouvierl (2003) study the production rate of BD s per star assuming 
BDs are a separate population, such as ejected embryos ( Reipurth &: Clarkd 
2001). Again they find that both, the physically very different environments of 
Taurus- Auriga and the ONC, can have produced the same ratios (about one BD 
per 4-5 stars) if BDs are ejected embryos with a dispersion of ejection velocities 
of about 2 km/s. 

Based on some additional observations, Luhman ( 2004! ) revised the Briceho et al 

(20021) results by finding that the number of BDs per star had been underes 



timated in Taurus- Auriga. Since the new spectroscopic study of ISlesnick et al 



(2004) also revised the number of BDs per star in the ONC downwards. iLuhman 



20041 ) retracts the significance of the claimed difference of the ratio in Taurus- 



Auriga and the ONC. Is a universal, invariant, BD production scenario still 
consistent with the updated numbers? 

Let the true ratio of the number of BDs per late-type star be 



R 



iV(0.02 -O.O8M ) _ JV] 
iV(0.15 - 1.0 M e 



BD.tot 



K 



(15) 



st, tot 



Note that here stars more massive than 10 M & are not counted because this 
would not make much sense given that Taurus- Auriga is mostly producing late- 
type stars given the limited gas mass available. But the observed ratio is 



R 



Ned. 



obs 



st,obs 



(1 + /) 



K 



st.tot 



(16) 



since the observed number of BDs, ABD,obs is the total number produced mul- 
tiplied by the fraction of BDs that are gravitationally bound to the popula- 
tion (B) plus the unbound fraction, U, which did not yet have enough time to 
leave the survey area. These fractions can be computed for dynamical mod- 
els of the Taurus-Auriga and ONC and depend on their mass and the dis- 
persion of velocities of the BDs. This velocity dispersion can either be the 
same as t hat of the stars if BDs form like stars, or larger if they are ejected 
embryos ( Reipurth &; Clarkd 200ll ). The observed number of "stars" is actu- 
ally the number of systems such that the total number of individual stars is 
Nst.tot = (1 + /) Ast i0bs , where / = N stjh - m /N st , sys is the binary fraction and 
^st.bin is the number of binary systems in a sample of iV stiSys stellar systems. 
Note th at here no distin ction is made between single or binary BDs. For Taurus- 
Auriga, Luhman ( 2004! ) observes i?TA,obs = 0.25 from which follows 



R TA = 0.18 since / TA = 1, B + U = 0.35 + 0.35 (17) 

( Kroupa et al J [20031 ). According to ISlesnick et al.l ( 2004! ) . the revised ratio for 
the ONC is i? NC obs = 0.28 so that 



# ONC = 0.19 because /qnc = 0.5, B + U = 1 + 



(18) 
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(|Kroupa et alJl2003h . Note that the regions around the stellar groupings in 



Taurus- Auriga not yet surveyed should contain about 30 % of all BDs, while all 
BDs are retained in the ONC. 

Therefore, the updated numbers imply that about one BD is produced per 
five late-type stars, and that the dispersion of eject ion velocities i s o~ P ,j ps 1.3 km/s. 
These numbers are an update of those given in iKroupa et ail (|2003l ). but the 
results have not changed much. Note that a BD with a mass of 0.06 M and 
a velocity of 1.3 km/s has a kinetic energy of 10~ 9 JW (pc/Myr) 2 which is 
rather comparable to the cut-off in BD-BD binding energies (Fig. [TH]) . This 
supports the notion that most BDs ma y be mildly ejected embry os. 

BDs can come in different flavours ( Kroupa &: Bouvierll2003l ) : star-like BDs, 
ejected embryos, collisional BDs and photo-evaporated BDs. As seen above, star- 
like BDs appear to be very rare, because BDs do not mix with stars in terms of 
pairing properties; some mechanism ensures the two types of object are kept sep- 
arated. Ejected embryos appear t o be a promising production channel, as noted 
above. IKroupa &: Bouvie r (2003) discount the collisional removal of accretion 
envelopes for the production of unfinished stars because this process is far too 
slow, or rare. The removal of accretion envelopes through photo-evaporation 
can occur, but o nly within the immediate vic inity of an O star and never 
i n Tau rus- Auriga ( Whitworth &: Zhmeckerj|2004l ). However. IKroupa &: Bouvier 
(2003) show that the radius within which photo-evaporation may be able to re- 
move substantial fractions of an accretion envelope within 10 5 yr is comparable 
to the cluster size in star-burst clusters that contain thousands of O stars. In 
such clusters photo-evaporated BDs may be very common. Globular clusters 
may then be full of BDs. 

Returning to the local star-forming regions, it thus now appears that the 
very different physical environments evident in Taurus- Auriga and the ONC pro- 
duce the same number of BDs per late-type star, so that there is no evidence for 
differences in the IMF across the hydrogen burning mass limit. There is also no 
substantial evidence for a difference i n the stell a r IMF in these two star-forming 
regions, contrary to the assertion by iLuhman (2004). Fig. \W\ shows the IMF 
in Taurus- Auriga: it is similar to the standard IMF (eq. [5D] below) , provided 
corrections for the high multiplicity are made. And in turn, the standard IMF 
is comparable to that seen in the ONC (Fig. [T6j) . The caveat that the classical 
pre-main sequence evolution tracks, upon which the observational mass deter- 
minations rely, are not really applicable for such young ages (§ 14. 2. p needs to be 
kept in mind though. 



6. The Shape of the IMF 

From the above discourse it thus becomes apparent that the distribution of 
stars by mass is a power-law with exponent or index a<i ps 2.35 for stellar masses 
m^O.5M . There exists a physical stellar ITlclSS limit, ?7imax* 

ps 150 M Q such 

that m < m max * (§ 13.1. p . The distribution of stars below the K/M dwarf tran- 
sition mass, 0.5 Mq, can also be described by a power law, but with exponent 
at ps 1.3 (§ S3) • Given the latest results described in § [571 the mass-distribution 
below the mass mi G (0.072 — 0.14) Mq is uncertain, but measurements indi- 
cate a power-law with exponent < ckq < 1- Because the binary properties of 
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Figure 19. The th ick histogram shows the Taurus- Auriga MF observed by 
iLuhman et al.l (|2003l ). and the dashe d curves are the stan dard IMF (eq. |20|) . 
The solid curves are the model from iKroupa et al.l (|2003f ) assuming all stars 
are in unresolved binary systems. The model shown by thick curves assumes 
stars only form in the mass interval 0.07 — 1.5 Mq, while the thin curves are 
for stars in the mass range 0.07 — 50 Mq. BDs would be part of the model as 
an extra population matched by the thick histogramme but are not explicitly 
shown in the model. The binary-star model (solid lines) is in g ood agreement 
with t he observed MF in the stellar regime. For details see IKroupa et al.l 
(2005). Note the discontinuity in the IMF evident between the dashed curve 
and the solid histogramme at 10 -12 Mq. 
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VLMSs and BDs differ substantially from those in the low-mass star regime, it 
can be argued that VLMSs and BDs need to be considered as a separate pop- 
ulation that is linked-to, but different from stars. This would then open the 
likely-hood that the IMF is discontinuous near mi, as may already be evident 
in Fig. [T9l And so fitting a functional description of the mass distribution with 
the continuity constraint would be wrong. For this reason it would appear to be 
dangerous, and in the worst case unphysical, to use one single function such as 
the log-normal form to attempt to describe the mass distribution from the frag- 
mentati on limit! 11 ! , mp ~ .01 M©, through to the physical stability limit m maxt . 
Indeed, IChabrien ( 2003a! ) already noted (i) that a single log-normal function 
over the entire stellar mass range is not compatible with the data, a power-law 
form for m > 1 Mq being required, and (ii) that the continuous IMF predicts 
too many VLMSs and BDs by a factor of about three. 

With these recent insights (power-law IMF over two orders of magnitude 
in mass and probable discontinuity near the sub-stellar mass limit), little of the 
argument for the advantages of a log-normal or any other mathematical form 
(Table[3]below) remains. Indeed, any such other mathematical form has the little 
nastiness that the tails of the distribution react to changes in the parametrisa- 
tion in a way perhaps not wanted when testing models. To give an example: 
a single log-normal form would change the slope of the IMF at large masses 
even if only the LF for late-type stars is to be varied. The standard (eq. I20|) 
two-part power-law stellar IMF, on the other hand, would allow changes to the 
index at low masses without affecting the high-mass end, and the addition of fur- 
ther power-law segments is mathematically convenient. The standard two-part 
power-law stellar IMF also captures the essence of the physics of star formation, 
namely a featureless power-law form for the largest range of stellar masses and 
a turnover near some fraction of a solar mass which probably is a result of in- 
creasing inefficiency of star formation (§173]). The overall shape of the IMF then 
sets the characteristic or average stellar mass. 



6.1. The standard or average IMF 

The various constraints arrived at above are summarised by a multi-part power- 
law, 



(mi) 



£(m) = k < 
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, mo < m < m± 
, mi < m < m,2 

, m n <m< m n+ i 



n = 0, 
71=1, 

n > 2, 
(19) 



When a cloud collapses its density increases but its temperature remains constant as long as 
the opacity remains low enough to enable the contraction work to be radiated away. The Jeans 
mass (eq. I14[l consequently decreases and further fragments with smaller masses form. When, 
however, the density increases to a level such that the cloud core becomes optically thick, then 
the temperature increases, and th e Jeans mass follows suit. Thus an opacity- limited minimum 
fragmentation mass is arrived at JLow fc Lvnden-BelllfilFFa iBoiail^ 
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where k contains the desired scaling and the mass-ratios ensure continuity, with 



a = +0.3 ±0.7 
ai = +1.3 ±0.5 
q 2 = +2.3 + 0.7 



0.01 < m/M Q < 0.08 
0.08 < m/M & < 0.50 
0.5 < m/M Q < 150 



n 
n 
n 



0, 

1, 
2. 



(20) 



The uncertainties are discussed in § 16.2.1 

This two-part power-law stellar IMF is referred to as the standard or canon- 
ical stellar IMF. It is also the average IMF since it rests upon a derivation from 
local star-counts for m < 1 M© that sample star-formation over the age of the 
MW disk (about 10 Gyr), and upon an averaging of the results obtained by the 
surveys of many OB associations and massive clusters for >fewM© (Fig. [1]). 
Eq. [20] is therefore a good bench-mark against which individually-derived IMFs 
can be compared to study possible systematic changes of the IMF with changing 
physical conditions. 

A larger value for ~ 2.7 in the stellar IMF may be closer to the truth than 
the above canonical S alpeter/Massey value ( o^ = 2.3) when binary-companions 
are corrected for (§[31 Sagar Richtlerlll99lh . in this case arrivina; at the KTG93 
IMF (IKroupa et alJll99£ 



q = +0.3 + 0.7 
ai = +1.3 + 0.5 
a 2 = +2.3 + 0.3 
a 3 = +2.7 + 0.7 



0.01 < m/M© < 0.08 , n = 0, 

0.08 < m/M© < 0.50 , n = 1, 

0.5 < m/M© < 1 , n = 2, 

1 < m/M Q , n = 3. 



(21) 



Note that in this case the observational data would not constrain the physical 
upper stellar mass limit ( Weidner fc Kroupal |2004| ). A detailed study of sys- 
tematic biases affecting power-law index measurements b ased on binned data is 



unde r way and will allow deeper insights into these issues (jMaiz Apellaniz et al 
20051 ). 

Note that 03 = 2.7 (the Scalo index) was also derived by IScalol ( 19861 ) for 



the MW disk population of massive stars such that the three-part-power law 
KTG93 IMF needs to be identified with the composite IMF introduced in § [87J 
Whether the stellar IMF needs to be changed away from the Salpeter/Massey 
value («3 = 2.3) because of corrections for unresolved multiples, and by how 
much the composite IMF then differs from the stellar IMF above 1M©, thus 
needs to be studied in more detail. 

6.2. The alpha plot 

A convenient way for summarising various studies of the IMF is by plotting 
indep endently-derived p ower-law in d ices in depende n ce of t he stellar mass ove r 
which they are fitted ([Scald Il998l : iKroupal l2001bl . 120021 : iHillenbrandl 120041 1. 
Fig. [20] plots such data: The shape of the IMF is mapped by plotting mea- 
surements of a at < Im >= {lm,2 — lm\)j2 obtained by fitting power-laws, 
£(m) oc m~ a , to logarithmic mass ranges lm\ to l m<i (not indic a ted here for clar - 
ity). Circles and triangles are data compiled by IScalol ( 19981 ); IKroupal ( 2001bl ) 
for MW and Large-Magellanic-Cloud (LMC) clusters and OB associations, as 
well as newer data, some of which are emphasised using different symbols (and 
colours) (Tabled]). Unresolved multiple systems are not corrected for in all these 
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data including the MW-bulge data. The standard stellar IMF (eq. |2"0"| corrected 
for unseen binary-star companions for m < 1M ) is the two-part power-law 
(thick short-dashed lines). 

Note that the M dwarf (0.1 — 0.5 Mq) MFs for the various clusters are 
systematically flatter (smaller ot\) than the standard IMF, which is mostly due to 
unresolved multiple systems in the ob served v a lues. S ome of the data do coincide 
with the standard IMF though, and iKroupal ( 2001bl ) argues that on correcting 



these for unresolved binaries the underlying true single-star IMF ought to have 
at ~ 1.8. This may indicate a systematic variation of ol\ with metallicity because 
the data are young clusters that are typically more metal-rich than the average 
Galactic field population for which ol\ = 1.3 (cf. eq. 1221 below) 

A power-law extension into the BD regime with a smaller index (ao = +0.3) 
is shown as a third thick short-dashed segment, but this part of the mass distri- 
bution may not be a continuous extension of the stellar distribution, as noted 
above. The upper and lower thin short-dashed lines are the estimated 99 % con- 
fidence range on on (eq. I20|) . Other binary-star-corrected solar-neighbourhood- 
IMF measurements are indicated as (magenta) dotted error-bars (Table [1]). Note 
that for m > 1 Mq correction for unseen companions may steepen the standard 
IMF toaRi 2.7. 

The long-dash-dotted horizontal lines in Fig. [20] labelled "SN" are those 
IMFs with as = 0.70(1.4) but Qo,ai,«2 as in eq. [201 for which 50 % of the 
stellar (including BD) mass is in stars with 8 — 50(8 — 120) Mq, respectively. It 
is noteworthy that none of the available clusters, not even including the star- 
burst clusters (Table [1]), has such a top-heavy IMF. Any hypothetical cluster 
so dominated by massive stars would disperse due to the mass loss from the 
supernova explosions. 

The vertical dotted lines in Fig. [20] delineate the four mass ranges (eq. [20]) , 
and the shaded areas highlight those stellar mass regions where the derivation of 
the IMF is additionally complicated especially for Galactic fi eld stars: for 0.08 < 
m/M (7) < 0.15 long pre- main sequence contraction times (jChabrier &: Baraffd 
20001 ) make the conversion from an empirical LF to an IMF (eq. [T]) dependent 



on the precise knowledge of stellar ages and the SFH, and for 0.8 < m/M Q < 2.5 
uncertain main-sequence evolution, Galactic-dis k age and the SFH of the MW 
disk do not allow accurate IMF determinations (Bi nnev et al. 20001 ). 



6.3. The distribution of data points in the alpha-plot 

The first thing to note about the data distribution in the alpha-plot is that 
there is (unfortunately) no clear systematic difference in IMF determinations 
with metallicity nor density of the population (cf. Fig. Q]). 

In order to understand the o rigin and nature of the dispersion of power-law 
indices evident in the alpha plot, IKroupal ( 2001bl ) investigates the dispersion of 



a values for a given mass interval theoretically. The result is that the dispersion 
can be u nderstood in terms of st atistical sampling from a universal IMF (as also 
found by lElmegreen 1997 . 19991 ) together with stellar-dynamical biases. Given 



such a theoretical investigation, it is now possible to compare the theoretical 
distribution of a values for an ensemble of star clusters with the observational 
data. This is done for stars with m > 2.5 Mq in Fig. [21] where the open (green) 
histogram shows the distribution of observational data from Fig. [20] The (blue) 
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Figure 20. The alpha plot (upper panel), and bolometric MLR, lL(lm), 
stellar main-sequence life-time, It, and possible range of Milky- Way (MW) 
disk ages shown as the shaded region (lower panel). Notation: Im = 
log lo (m/M ), It = log 10 (r/yr), IL = log 10 (L/Lo). For more details see 
text. 



48 



Pavel Kroupa 



shaded histogram is the theoretical ensemble of 12 star clusters containing ini- 
tially 800 to 10 4 stars that are "observed" at 3 and 70 Myr: stellar companions 
are merged to give the system MFs, which are used to measure a, but the input 
single-star IMF is in all cases the standard form (eq. I20|) . The dotted curves are 
Gaussians with the same average a and standard deviation in a obtained from 
the histograms. Fixing at =<a> and using only| a\< 2a a for the observational 
data gives the narrow thin (red) dotted Gaussian distribution which describes 
the Salpeter peak extremely well (not by construction). 

The interesting finding is thus that the observational data have a very pro- 
nounced Salpeter/Massey peak, with broad near-symmetric wings. This indi- 
cates that there are no serious biases that should skew the distribution. For 
example, if the observational sample contained clusters or associations within 
which the OB stars have a low binary fraction compared to others that have a 
very high multiplicity fraction, we would expect the binary-deficient cases to de- 
viate towards low a values. Likewise, energetic dynamical ejections from cluster 
cores would skew the IMF towards more massive stars (reduced a) while mass 
segregation has a similar effect. 

In contrast, the theoretical data show (i) a distribution with a mean shifted 
to smaller a ~ 2.2 that has (ii) a larger width than the observational one. The in- 
put standard Salpeter/Massey index is not really evident in the theoretical data, 
and if these were the observational data then it is likely that the astronomical 
community would strongly argue for the case that the IMF shows appreciable 
variations. It is peculiar that the empirical data are so much better behaved, 
since all the additional complications (theoretical stellar models, rotating/non- 
rotating stars) are not modelled in the theoretical data. The elucidation of the 
reason for the difference between the much more "well-behaved" observational 
data and the theoretical data will need further theoretical work which will have 
to attempt to re-produce the observational procedure as exactly as is possible. 



6.4. Mass density and some other numbers 

The nearby Hipparcos LF, ^near 

(Hipp), has p = (5.9 ± 0.3) x 10 3 stars/pc 3 in 
the interv al My = 5.5 — 7.5 corresponding to the mass interval ni2 = 0.891 — 
0.687 M (|Kroupal [20013 ^ using the KTG93 MLR (Fig. EJ. f™*£(m)dm = p 

yields k = 0.877 ± 0.045 stars/ (pc 3 M ). The number fractions, mass fractions 
and Galactic-field mass densities contributed by stars in different mass-ranges 
are summarised in Table [2j 

The local mass density made up of interstellar m atter, p gas rs 0.0 4 ± 
O.O2M /pc 3 , and stell ar remnants, p re m » 0.003 M /pc 3 (lWeidemannlll990l or 



pvem ^ 0.005 M /pc 3 ( Chabri er 2003a | and references th erein) . Giant stars con- 



tribute about 0.6 x 10~ 3 M pc~ 3 (jHavwood et al.lll997l ). so that main-sequence 
stars make up about half of the baryonic matter density in the local Galac- 
tic disk (Table [2]). BDs, which for some time were regarded as candidates for 
contributing to the dark-matter problem, do not constitute a dynamically im- 
portant mass component of the Galaxy, even when eq. [20] is extrapolated to 
0.0 Mq giving p BD = 3.3 x lO~ 3 M /pc 3 . This is corroborated by dynamical 
analysis of local stellar space motions that imply there is no dark matter in the 
Milky- Way disk dFlvnn k Fuchslll994l ). and the revision of the thick-disk mass 
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Figure 21. Distribution of a values for massive stars. See text for details. 



50 



Pavel Kroupa 



density to larger values ( Fuhrmann 2004 ; Soubiran et al. 2003t) fur ther reduces 



the need for dark matter within the solar circle (cf IChabrierl 2003al ) . 



Table [2] also shows that a star cluster looses about 13 per cent of its mass 
through stellar evolution within 10 Myr if 03 = 2.3 (turnoff-mass mto ~ 20 Mq), 
or within 300 Myr if 03 = 2.7 (turnoff-mass m to ~ 3M@). After 5 Gyr the mass 
loss through stellar evolution alone amounts to about 45 per cent if «3 = 2.3 or 
29 per cent if 03 = 2.7. Mass loss through stellar evolution therefore poses no 
risk for the survival of star clusters for the IMFs discussed here, since the mass- 
loss rate is small enough for the cluster to adiabatically re-adjust. A star-cluster 
would be threatened through mass loss from supernova explosions if q ~ 1.4 for 
8 < m/ Mq < 120 which would mean a mass- loss of 50 per cent within about 
40 Myr when the last supernova explodes. It is remarkable that none of the 
measurements has found such a low a for massive stars (Fig. [20]) . 

6.5. Other IMF forms 

As already discussed at the beginning of this section, the standard or canonical 
power-law IMF (eq. [20]) provides a good description of the data combined with 
mathematical ease and physical meaning. A strong advantage of this parametri- 
sation is that each section can be changed without affecting another part of the 
IMF, as stated above. As an explicit example, should the BD MF be revised, ao 
can be adopted accordingly without affecting the rest of the mass distribution 
in the well-tested stellar regime. 

But additional forms are in use and are preferred for some investigations: 
In Fig. [20] the quasi-diagonal (black) lines are alternative analytical forms sum- 
marised in Table [3] Of these the generalised Rosin- Rammler function (eq. Ch in 
Table [3] thick short-dash-dotte d curve) best rep resents the data, apart from a 
deviation for m<^ 10 Mq, which IChabrierl ( 2003al ) fixes by adopting a Salpeter/- 



Massey power-law extension for m > IMq. Interpretation of m in terms of a 
characteristic mass poses difficulties. The popular Miller-Scalo log-normal IMF 
(eq. MS in Table [3]) deviates even more strongly from the empirical data at 
high masses. Larson's eq. Lb in Table [3] fits rather well, except that it may 
predict too many BDs. Finally, th e effec tive initial mass function for galactic 



disks proposed by iHollenbach et al.l (120051) feci. Holt) r eproduces the data in the 



alpha-plot quite well (fig. 1 in Hollenbach et al. 20051 ). and is not incorporated 



into Fig. [20] here. Note however that a composite IMF ought to be steeper (have 
a larger a) at high masses (§E]). 

The closed functional IMF formulations (eqs. MS, La, Lb, Ch, Holt) have 
the advantage that possible variations of the IMF with physical conditions can be 
studied more naturally than with a multi-power-law form, because they typically 
have a characteristic mass that can be varied directly. 



7. Variation of the IMF and theoretical aspects 

Is the scatter of data points in the alpha-plot (Fig. l20|) a result of IMF variations? 
For this to be conclusively convincing would require a measurement to lie further 
from the standard IMF than the conservative uncertainty range adopted in eq. [20l 
and shown in the figure. However, the adopted uncertainties on in eq. 1201 stem 
from the scatter in the alpha plot, so that this argument is circular. 
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An independent indication of the uncertainties inherent to IMF determi- 
nations can be obtained by comparing IMF estimates of the same cluster by 
different authors. This is demonstrated for the well-studied Pleiades, ONC and 
for 30 Dor (Table [JJ Fig. [20]) . Overall, the uncertainties in a are about ±0.5 
which is also about the scatter evident in all the data, so that there is no indica- 
tion of significant outliers (except in the shaded regions, see below). Substantial 
differences for VLMSs and BDs are evident for the extremely young ONC al- 
lowing an estimate of likely non-physical variations in the alpha plot. Data 
reduction at these l ow masses is hamp ered by variability, variable reddening, 
spurious detections (jSlesnick et al.|[2004 ). It is clear that because the procedure 
of measuring a(lm) is not standardised and because the IMF is not a single 
power-law, author-author variations occur simply due t o the use of different 
mass ranges when fitting power-laws. Here the work by iMaiz Apellaniz et al. 
(200§|), who is addressing such issues, is important. 

Significant departures from the standard IMF only occur in the shaded ar- 
eas of the alpha plot. These are, however, not reliable. The upper mass range 
in the shaded area near 1 Mq poses the problem that the star-clusters inves- 
tigated have evolved such that the turn-off mass is near to this range. Some 
clusters such as p Oph are so sparse that more massive stars did not form. In 
these cases the shaded range is close to the upper mass limit leading to possi- 
ble stochastic stellar-dynamical biases since the most massive stars meet near 
the core of a cluster due to mass segregation, but three-body or higher-order 
encounters there can cause expulsions from the cluster. Furthermore, p Oph is 
still forming, leading to unknown effects that are likely to enhance variations 
in the first derivative of the IMF (i.e. in a values). Dynamical ejections are 
probably the cause for the "OB field-star" MF which has a ~ 4.5 and is in- 
terpreted to b e the result of isolated high-mass star-formation in small clouds 
(|Massevll 19981 ). Precise proper-motion measurements have, however, shown that 
even t he best candidates for this isolated population have very high space mo- 
tions ( Ramspeck et al1l200ll ) which are the result of energetic stellar-dynamical 
ejections when massive binary systems interact in the cores of star-clusters in 
normal but intense star-forming region s, thus posing important constraints on 
the properties of OB binary syste ms ( Clarke <fe Pringief 19921 ; Kroupa 2001a : 



Pflamm-Altenburg &: Kroupall2006l ). 



The shaded area near 0.1 Mq poses the problem that the VLMSs are not on 
the main sequence for most of the clusters studied, and are again prone to bias 
through mass-segregation by being underrepresented within the central cluster 
area that is easiest to study observationally. Especially the latter is probably 
biasing the M35 datum, but some effect with metallicity may be operating, 
especially so since M35 appears to have a smaller a near the H-burning mass 
limit than the Pleiades cluster which has a similar age but has a larger abundance 
of metals (Fig. [To]) . The M35 cluster ought to be looked at again with telescopes. 

To address stellar-dynamical bia ses, an extensiv e theoretical library of binary- 
rich star clusters has been assembled ( Kroupall2001bl ) covering 150 Myr of stellar- 
dynamical evolution taking into account stellar evolution and assumin g the stan- 
dard IMF (eq. [20]) in all cases. This extends the notion raised by lElmegre"en 
(1999) that the scatter seen in the empirical alpha-plot may be due to finite— N 
sampling from a universal IMF that can be interpreted as a probability density 
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distribution. Evaluating the MF within and outside of the clusters, at different 
times and for clusters containing initially 800 — 10 4 stars, leads to a theoretical 
alpha-plot which reproduces the spread in a(lm) values evident in the empir- 
ical alpha-plot quite well. This verifies the conservative uncertainties adopted 
in the standard IMF (eq. [20]) but also implies that the scatter in the empirical 
alpha-plot around the standard IMF cannot be interpreted as true variations. 
This work also addresses the bias due to unresolved binary systems, with the 
result that typical open clusters ought to be underrepresented in VLMSs. This 
bias is incorporated in the definition of the standard IMF, as is evident in the 
empirical alpha-plot. 

One peculiar feature of the empirical alpha-plot that merits further study 
is the distribution of a data for m > 2.5 M & (§ 16.3. p . A histogram of the 
empirical data (Fig. [2T]) shows a narrow peak positioned at the Salpeter value, 
with symmetric broad wings. The empirical data are not distributed normally. 
The aforementioned theoretical alpha-plot shows a very different distribution, 
the width being larger than for the empirical data. Interestingly, the spread, 
a a,f = 0.08, of the narrow peak in the empirical data is very similar to the 
uncertainties quoted by Massey in an extensive and consistent observational 
determination of the IMF for massive stars, a = 2.2±0.1 (Fig.[T]). If a = 2.3±0.1 
is adopted for massive stars, then the measurement a = 1.6±0.1 for the massive 
but difficult-to-observe Arches cluster near the Galactic centre (Table [T]) would 
definitely mean an IMF that is top-heavy for this extreme population. But as 
discussed in § [Lj the Arches is suffering heavily from the tidal forces that strip its 
less massiv e mem bers so that the PDMF is likely skewed towards massive stars. 
Kim et al. (|2006l ) indeed verify the near-Salpeter/Massey IMF in the Arches 



cluster. The question thus remains if the non- Gaussian distribution of empirical 
a values contains information on a possible physical variation of the IMF. 

Two other well-studied massive star-burst clusters have a ~ Salpeter /Massey 
(30 Dor and NGC 3603, Table [1]) implying no clear evidence for a bias that re- 
solved star-burst clusters prefer smaller a and thus more massive stars relatively 
to the number of low -mass stars. Low-mass stars are known to form in 30 Dor 
( Sirianni et al. 2000l ). although their MF has not been measured yet due to the 
large distance of about 55 kpc. From the ONC we know that the entire mass 
spectru m 0.05 ^m/M^ ^60 is represented r oughly following the standard IMF 
(Fig. 116] Pflamm-Altenburg &: Kroupa 20061 ). The Pleiades appears to have had 
an IMF very similar to the standard one (Fig. fT5"|) . although for massive stars 
a steeper IMF with a 3 ~ 2.7 appears to be suggested by the theoretical work 
(jMoraux et al]l2004T l. 



The available evidence is thus that low-mass stars and massive stars form 
together even in extreme environments without, as yet, convincing demonstra- 
tion of a variation of the number ratio. This is also supp orted by an impressive 
observational study ( Luhman et al. 2000l ; Luhman 20041 ) of many close- by star- 
forming regions using one consistent methodology to avoid author-author varia- 
tions. The result is that the IMF does not show measurable differences from low- 
density star-forming regions in small molecular clouds (n = 0.2 — 1 stars/pc 3 in 
p Oph) to high-density cases in giant molecular clouds (n = (1—5) x 10 4 stars/pc 3 
in the ONC). This result extends to the populations in the truly exotic ancient 
and metal-poor dwarf-spheroidal satellite galaxies which are speculated to be 
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severely dominated by dark matter but definitely constitute star-forming con- 
ditions very different from present-d ay events. Two such close companions to 
the Milky Way have been observed ( Grillmair et al.|[l998l ; iFeltzing et al.l fl999l ) 
finding the same MF as in globular clusters for 0.5 ^m/M ^0.9 and thus no 
very evident differences to the standard IMF. 

However, there are peculiar indications of top-heavy IMFs such as in some 
massive star-burst clusters in the M82 galaxy. Spectroscopy of the unresolved 
M82-F cluster derives, via the inferred velocity dispersion, a mass and from the 
luminosity a mass-to-light ratio that is significantly smaller than the ratio ex- 
pected from the standard IMF for such a young (about 60 Myr) population. 
The implication is that the M82-F population may be significantly depleted in 
low-mass stars, provided the velocity dispersion is representative of the entire 
cluster. A possibility that will have to be addressed using stellar-dynamical 
modelling of forming star clusters is that M82- F may have lost low-mass stars 
due to tidal shocking ( Smith &; Gallagher] l200ll ). Highly pronounced mass seg- 
regation which leads to a dynamically decoupled central core of O B stars is an 



impo rtant mechanism for reducing the measured mass-to- light ratio (jBoilv et al 



2005), while rapid expulsion o f residual gas from forming clusters enhances the 



measured mass-to-light r atios (jGoodwin <k Bastianll2006l ). 

The proclamation bv lBriceho et al] ( 20021 ) that Taurus- Auriga is producing 
significantly fewer BDs per star than the ONC, and thus for a possible first robust 
indication of a varia b le IMF, was found to not be the correct interpre tation by 
iKroupa et al.1 (|2003h : IKroupa fc Bouvierl <|2003T > and was retracted by iLuhmar] 
(2004) based on newer data. A new analysis of these actually shows that Taurus- 
Auriga and the ONC are producing the same number of BDs per star, about 0.2 
(§ Nevertheless, the discussion in § ED shows that BDs probably need to 
be treated as a separate population such that the IMF may be discontinuous 
near the stellar /sub-stell ar mass limit. Since BDs may come in four flavours, 
Kroupa Sz Bouvierl ( 20031 ) suggest that the photo-evaporated sort could far out- 



weigh the others in star-burst clusters, so that globular clusters could contain a 
larger fraction of BDs per star than the current modest Galactic clusters. 

Differences in the metallicity, Z, of the popul ation also do not lead to ob- 
servable variations of the IMF for massive stars (|Massevi I2003L Fig. HJ. The 
mass of the most massive star, m max * « 150 M , is independent of Z. Thus 
the distribution of masses for massive stars does not appear to be affected by 
the metallicity of the star-forming gas and therefore radiation pressure on dust 
grains during star-assembly cannot be a physical mechanism establishing m max * 

(§©■ 

However, there may be an effect for VLMSs. Present-day star-forming 
clouds typically have somewhat higher metal-abundances ( [ Fe/H l~ +0.2) com- 
pared to 5 Gyr ago ([Fe/H]« -0.3) (|Binnev Merrifieldl fl998l ) which is the 
mean age of the population defining the standard IMF. The data in the empiri- 
cal alpha-plot indicate that some of the younger clusters may have a single-star 
MF that is somewhat steeper (larger a±) than the standard IMF if unresolved 
binary-stars are corrected for (Fig. [20]) . This may mean that clou ds with a 
larger [Fe/H] produce relatively more low- mass stars ( Kroupa 2001bl ) which is 
tentatively su pported by the typicall y but not significantly flatter MFs in glob- 
ular clusters (jPiotto Zocca li 1999) that have [Fe/H] : 



-1.5, and the recent 
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finding that the old and metal-poor ([Fe/H]= —0.6) thick-disk popu lation also 
has a flatter MF below 0.3 M with ot\ « 0.5 (|Revle fc Robinll200ll ). If such a 
systematic effect is present, then for m£O.7M and to first order, 



am 1.3 + Aq [Fe/H], 



(22) 



with Aa ~ 0.5, being similar to the adopted uncertainty on a. Theoretical 
considerations do suggest that for sufficiently small metallicity a gas cloud cannot 
cool efficiently causing the Jeans mass required for gravitational collapse to 
be larger. In particular, the first stars ought to have large masses owing to 
this effe ct and the generally higher ambie nt temperatures at early cosmological 
epochs (|Larsonlll998l : H romm et al.ll200ll ). Finding the remnants of these poses 
a major challenge. 

An easier target is measuring the IMF for low-mass and VLMSs in metal- 
poor environments, such as young star-clusters in the Small Magellanic Cloud. 
That metallicity does play a role is becoming increasingly evident in the planetary- 
mass regime in that the detecte d exo-planets occur mostly around stars that are 
more metal-rich than the Sun (jVogt et al.ll2002l ; IZucker L Mazehll200ll ). 

While the Jeans-mass argument (eq. [T4"|) should be useful as a general indi- 
cation of the rough mass scale where fragmentation of a contracting gas cloud 
occurs, the concept clearly breaks down when considering the stellar masses 
that form in star clusters. The central regions of these are denser, formally 
leading to smaller Jeans masses which is the opposite of the observed trend, 
where even in very young clusters massive stars tend to be located in the inner 
regions. More complex physics is clearly in volved. Stars may, to a c ertain ex- 
tend, regulate their own mass by outflows ( Adams &: Laughlin 1996), and the 
coagulation of proto-stars probably plays a significant role in the densest regions 
where the cloud-core collapse time, r co ii, is longer than the fragment collision 
time-scale, r cr . The collapse of a fragment to a proto-star takes no longer than 
r co ii « 0.1 Myr (|Wuchterl fc Klessenl l200"lh . so that 



WMyr ~ 42 



(M ecl /M ; 



< 0.1 Myr, 



(23) 



where R is the half-mass radius of a Plummer-sphere model, implies M/R 3 > 
lO 5 M0/pc -3 . Such densities, where proto-stellar interactions are expected to 
affect the emerging stellar mass-spectrum, are found in the centres of very dense 
and rich embedded s tar clusters such as the ONC before they expand as a re- 
sult of gas expulsion ( Kroupa et al.|[200ll ; IVine Bonnelll 12003). Thus, only for 
massive stars is the form of the IMF probably affected by coagulation, which 
may expl ain why massive stars are usually centrally concentrated in very young 
clusters ( Bonnell et al,lfl998l ; lKlessenll200ll ). However, the observed mass seg- 
regation in very young clusters cannot as yet be taken as conclusive evidence 
for primordial mass segregation and coagulation, unless precise iV-body compu- 
tations of the embedded cluster are performed for each case in question. For 
example, models of the ONC show that the degree of observed mass segregation 
can be established dynamically within about 2 Myr (Fig. [T7|) despite the em- 
bedded and much denser configuration having no initial mass segregation. The 
notion be hind such an assumption is that star clusters fragment h eavily sub- 
clustered ( Megeath et al.|[l996l ; iBontemps et all l2Q0ll : iKlessenlbOQll ). and each 
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sub-cluster may form a few OB stars with a few hundred associated lower-mass 
stars (Table [2]) , so that the overall morphology may resemble a system without 
significant initial mass segregation. The theoretical time-scale, ti — 1\ in Fig. [T71 
for mass segregation to be established can be shortened by decreasing the relax- 
ation time. This can be achieved by reducing the number of stars in the model, 
for example. But it may prove impossible to find agreement at the same time 
with the density profile and kinematics because the ONC is probably expanding 
rapidly now. Clearly this issue requires more study. 

The origin of most stellar masses is indicated by a remarkable discovery for 
the low-mass p Oph cluster in which star-formation is on-going. Here the pre- 
stellar and proto-star MF are indistinguishable and both are startlingly similar 
to the st andard IMF, even i n showing the same fl attening of the power-law at 
0.55 Mq dMotte et al.lll998l : iBontemps et al.1l200lh . The pre-stellar cores have 
sizes and densities that are in agreement with the Jeans-instability argument 
for the conditions in the p Oph cloud, so that cloud-fragmentation due to the 
collapse of density flu ctua tions in a dissipating turb ulent interstellar medium 
dNordlund fc Padoan|[2003l:lPadoan fc Nordlundll2002l :lMac Low fc Klessenll2004l : 



Tillev Pudrital2004l ; lPadoan &: Nordlundll2004l ) appears to be the most-important 



mechanism shaping the stellar IMF for masses 0.5 ^m/MQ^S a few Mq, the 
shape of the IMF being determined by the spectrum of density fluctuations 
in the molecular cloud. Similar results h ave been obtained for the Serpens 
cloud s and for the clouds in Taurus-Auriga ( Testi fe Sargentlll998l ; lOnishi et al 



20021 . respectively). The majority of stellar masses making up the standard 



IMF thus do not appe ar to suffer significa nt subsequent modifications such as 
competitive accretion ( Bonnell et al. 200 ll ). proto- stellar mergers or even self - 
limitation through feedback processes. The work of iPadoan fc Nordhmdl ([2002) 
demonstrates that, under certain reasonable assumptions, the mass function 
of gravitationally unstable cloud cores deriving from the power-spectrum of a 
super-sonic turbulent medium leads to the observed standard IMF above IMq. 
The flattening at lower masses is a result of a reduction of the star-formation 
efficiency because at small masses only the densest cores can survive sufficiently 
long to collapse. 

The intriguing result from p Oph, Serpens and Taurus-Auriga, in which the 
stellar, proto-stellar and pre-stellar clump mass spectra are similar to the stellar 
IMF (eq. [20]) . is consistent with the independent finding that the properties of 
binary systems in the Galactic field can be understood if most stars formed in 
modest, p Oph-t ype clusters w ith primordial binary properties as observed in 
Taurus- Auriga ( Kroupalll995bl ). and with the independent result derived from 
an analysis of the distribution of local star clusters th at most stars appear to 
stem from such modest clusters ( Adams fc Mversll200ll ). However, the standard 
IMF is also similar to the IMF in the ONC (Fig. I16j) implying that fragmentation 
of the pre-cluster cloud proceeded si milarly. 

The impressive computations by I Bonnell Sz Bate] ( 20021 ) and collaborators of 
dense clusters indeed not only predict the observed m max (M ec i) relation (Fig. [5]), 
but they also show that a Salpeter/Massey power-law IMF is obtained as a result 
of competitive accretion and the merging of stars near the cluster core driven by 
accretion onto it. The reason as to why the IMF is so invariant above a few M & 
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may thus be that the various physical processes all conspire to give the same 
overall scale- free result. 

This still leaves the origin and nature of VLMSs and BDs unclear. As the 
discussion in § [5] suggests, VLMSs and BDs probably need to be considered 
a different or extra population whic h does not mix well, in terms of pairing, 



with stars. iMaxted &; Jeffries! (|2005l ) demonstrate that virtually all currently 



available theoretical work on the formation of BDs is excluded by the data. To- 
gether with the realisation that existi ng star-formation theory fails to reproduce 



the b inary properties of young stars ([Goodwin et al.ll2004l ; iGoodwin h Kroupa 



20051 ) . this suggests that our theories lack major ingredients that probably are 



related to stellar feedback processes. That VLMS and BD binaries have a sig- 
nificantly higher energy scale than stellar binaries does, however, suggest that 
their formation may be linked to a dense environment which they probably leave 
abruptly. 



8. Composite stellar population 

We have thus seen that while a conclusive theoretical formulation of the IMF is 
still wanting, we do have some good ideas about its origin and a good impression 
of its shape, and empirically there is not a strong case for systematic variation 
of the IMF with physical conditions. It is thus reasonable to assume that the 
stellar IMF is invariant. 

The natural assumption has often been made that independent of the star- 
formation rnodej_Jl^steJIar distribution is sampled randomly from the invariant 
IMF (e.g. Elmegreen 20041 ). Thus, for example, 10 5 clusters, each with 20 stars, 



would have the same composite (i.e. combined) IMF as one cluster with 2 x 
10 6 stars. 

However, the existence of the m max (M ec i) relation (§ 13.1. p has profound 
consequences for composite populations. It immediately implies, for example, 
that 10 5 clusters, each with 20 stars, cannot have the same composite (i.e. com- 
bined) IMF as one cluster with 2 x 10 6 stars, because the small clusters can never 
make stars more massive than about 1.5 AT© (Fig. [5]). Thus, galaxies, that are 
composite stellar populations consisting of many star clusters, most of which 
may be dissolved, would have steeper composite, or integra ted galactic IMFs 
(IGIMFs), than the stel lar IMF in each individual cluster ( Vanbeverenl Il982l ; 
Kroupa fc Weidneril2003n . 



The IGIMF is an integral over all star-formation events in a given star- 
formation "epoch" t,t + 5t, 

(SFR(t)) 

f igimf ("»; t)= £ (m < m max (M ed )) £ cc i(M ccl ) dM ecl . (24) 

J Mad ,min 



Here £(m < m max ) £, ec i(M ec \) dM cc \ is the stellar IMF contributed by £ ec i dM ec i oc 
MZa dM ec i clusters with mass near M ec \. M ec i jmax follows from the maximum 
star-cluster-mass vs global-star-formation-rate-of-the-galaxy relation, 

M ecl , max = fn(SFR) (25) 
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(eq. 1 in lWeidner & Kroupall2005l . as derived bv lWeidner et al.ll2003 ). M ecl min 



5 Mq is adopted in the standard modelling and corresponds to the smallest star- 
cluster units observed. At time t the SFR is 

SFR(t) = (26) 
ot 



where 



M tot = f CCl ' maX M ecl £ ecl (M ecl )dM ecl (27) 



is the total stellar mass assembled in time 5t which Weidner et al.l ( 20041 ) de- 



fine to be a "star-formation epoch", within which the ECMF is sampled to 
completion. This formulation leads naturally to the observed M ec i jmaiX (SFR) 
correlation if the ECMF is invariant, (3 ~ 2.35 and if the "epoch" lasts about 
5t = 10 Myr. Thus, the embedded cluster mass function is fully sampled in 
10 Myr intervals, independent of the SFR. This time-scale compares very well 
indeed to the star -formation time-scale in normal galactic disks measured by 
Egusa et al. ( 20041 ) using an entirely independent method, namely from the 



offset of HII regions from the molecular clouds in spiral- wave patterns. The 
time-integrated IGIMF then follows from 

£igimf(«"i) = / Zigimf (rn;t)dt, (28) 
J o 

where tq is the age of the galaxy under scrutiny. 

Note that £iGiMF( m ) is the mass function of all stars ever to have formed 
in a galaxy, and can be used to estimate the total number of supernovae ever to 
have occurred, for example. £lGlMF(m; t), on the other hand, includes the time- 
dependence through a dependency on SFR(t) of a galaxy and allows one to 
compute the time-dependent evolution of a stellar population over the life-time 
of a galaxy. 

Because more-massive stellar clusters are observed to form for higher SFRs 
(eq. [25]) . the ECMF is sampled to larger masses in galaxies that are experiencing 
high SFRs, leading to IGIMFs that are flatter than for low-mass galaxi e s tha t 
have had only a low-level of star-formation activity. Weidner fc Kroupa ( 20051 ) 



show that the sensitivity of the IGIMF power-law index for m ^ 1 M & towards 
SFR variations increases with decreasing SFR. 

Thus, galaxies with a small mass in stars can either form with a very 
low continuous SFR (appearing today as low-surface-brightness but gas-rich 
galaxies) or with a brief initial SF burst (dE or dSph galaxies). The IGIMF 
ought to vary significantly among such galaxies (Fig. [22]) . In all cases, however, 
the IGIMFs are invariant for m^l.3M Q which is the maximal stellar mass 
in 5Mq "clusters" (Fig. [5]). Low-surface-brightness galaxies would therefore 
appear chemically young, while the dispersion in che mical properties ought to 



be larger for dwarf galaxies than for larger galaxies ([Goodwin Page 



>u gnt to 
si 120051 : 



Weidner &: Kroupall2005l ; iKoppen et al] [2007). Another interesting implication 



is that the number of supernovae per star would be significantl y smaller over cos- 
mological times than predicted by an invariant Salpeter IMF ( Goodwin &: Page! 
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Figure 22. The power-law index ckigimf of the IGIMF above 1.3 Mq as a 
function of the stellar mass of a galaxy, which determines the average SFR 
over a Hubble time. The lower bound is for an initial SF burst that forms 
the entire stellar galaxy, while the upper bound is for a constant SFR over a 
Hubble time. In all cases the input stellar IMF is the standard form (eq. |20|) . 
Note that the detailed results d epend on the adopted form of the ECMF (cf. 
Fig. [231). For further details see lWeidner fc Kroupal (|2005| ). 
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Figure 23. The number of supernovae per star in the IGIMF divided by 
the number of supernovae per star in the standard IMF, 77, as a function of 
the stellar galaxy mass. The upper shaded area is for an ECMF with (3 = 2. 
while the lower shaded area assumes = 2.35. The upper bound for each 
shaded region is for an initial SF burst model, while the lo wer bounds are for a 
consta nt SFR over a Hubble time. For further details see lWeidner &: Kroupal 

IHl). 
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These new insights should lead to a revision of theoretical work on galaxy 
formation that typically until now relied on an invariant IMF. Empirical evidence 
in favour of or against the notion of a galaxy-variable IGIMF is being studied 
and will ultimately lead to a refinement of the ideas. 

9. Concluding comments 

Spectacular advances have been achieved over the past decade in the field of 
IMF research and this affects a vast area of astrophysics. 

The stellar IMF appears to be extra-ordinarily universal. It is described by 
eq. [20j The evidence for top-heavy IMFs comes from either unresolved clusters 
or clusters that are very difficult to observe. The most significant uncertainty 
that remains at high masses despite the presence of exquisite data is the exact 
shape of the IMF for massive stars, because the biases due to unresolved multiple 
systems and due to stellar rotation have not been studied in sufficient detail. 
The true IMF may be closer to the Scalo value «3 ~ 2.7 rather than being 
Salpeter/Massey. For sub-stellar masses more data are necessary to improve 
constraints on ao and to better quantify any claims for a variation of the IMF. 
Further problems are inadequate stellar models for stars younger than a few Myr, 
and the under-representation of stellar-dynamical modelling which, however, is 
an absolute necessity in the search for IMF variations among clusters. 

Among intriguing recent results are that the BDs appear to be a distinct 
population from that of low-mass stars; their pairing properties have a different 
energy scale. This sets the stage for a probable discontinuity in the IMF near 
the stellar /sub-stellar boundary which has already probably been detected. This 
then questions the validity of fitting continuous log-normal functions across the 
VLMS/BD barrier. Furthermore, the IMF does appear to have a physical max- 
imum stellar mass that has now been found empirically. Stars with m ^ 150 M@ 
do not appear to exist, unless they implode invisibly shortly after being formed. 
In addition, by realising that star clusters are the true fundamental building 
blocks of a galaxy, we merely need to add up all clusters and their (invariant) 
IMFs to arrive at the integrated galactic initial mass function. The IGIMF 
varies according to the star-formation history of the galaxy. This formulation 
now allows us to compute the IGIMF, as a function of time, for galaxies that are 
converting their gas supply only slowly to stars, or for galaxies that like to eat 
it up all at once. The implications are that equally-old galaxies can have very 
different chemical compositions ranging from young to evolved, and that the 
cosmological supernova II rate may be significantly smaller and dependent on 
galaxy type than if an invariant Salpeter IMF is assumed. Galaxy-formation and 
evolution computations with this latter assumption are not likely to be correct. 

Thus, the IMF research as presented here has advanced onto cosmological is- 
sues, whereby many details still need to be worked out. But the most intimately 
connected research field, namely star formation, is being tested rather dramat- 
ically. The brilliant computational results that have been becoming available 
have lead to a beautiful agreement with the observed maximum-stellar-mass vs 
cluster-mass correlation and reproduce the general shape of the IMF. However, 
different ideas (competitive accretion, coagulation and simply the distribution 
of gravitationally unstable regions in turbulent clouds) all lead to virtually the 
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same type of theoretical IMF. Detailed computations of the formation of stars 
actually get it all wrong - the primordial-binary properties cannot yet be de- 
scribed, and the failure of modern theory is especially evident in the BD regime, 
where virtually all current ideas are being excluded by the data. The failures of 
modern theory are not surprising though, because the primary process driving 
star-formation is injection of feedback energy, the consistent inclusion of which 
goes beyond currently available computational resources. On the bright side, the 
convergence of different mechanisms to the same Salpeter/Massey IMF may be 
interpreted to simply mean that the shape of the IMF for m ^ 0.5 M is indeed 
invariant to the physics, as is observed. 
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a 


a 




mass range [M©] 


mass range [M©] 


mass range [M©] 


Orion Nebula cluster, ONC 








Muench et al. (2000) 


-0.35 


+ 1.25 


+2.35 


magenta small open circles with central dot 


0.02-0.08 


0.08 - 0.80 


0.80 - 63.1 


magenta large open circles with central dot 


+0.00 
0.02-0.08 


+ 1.00 
0.08 - 0.40 


+2.00 
0.4-63.10 


Hillenbrand & Carpenter (2000) (HC00) 
magenta large thick open circle 
with central dot 


+0.43 
0.02 - 0.15 






Luhman et al. (2000) 


+0.70 






magenta small thick open circle 
with central dot 


0.035 - 0.56 






Pleiades 








Moraux et al. ('2001) 


+0.51 ± 0.15 






green circles with central dot 


0.04 - 0.30 






Hamblv et al. (1999). from Barrado v Navascues et al. (2001) 


+0.56 


+2.67 




green circles with central dot 


0.065 - 0.60 


0.6 - 10.0 




a Ori 








Beiar et al. (2001) 


0.8 + 0.4 






green solid circle 


0.013 - 0.20 






M35 








Barrado v Navascues et al. (2001) 


-0.88 + 0.12 


0.81 + 0.02 


2.59 + 0.04 


T ■ 7 7 1 

green solid circle 


0.08 - 0.2 


0.2 - 0.8 


0.8 - 6.0 


IC 348 








Naiita et al. (2000) for MLR from Baraffe et al. (1998) 


+0.5 






green solid circle 


0.015-0.22 






NGC 2264 








Park et al. (2000) 






+2.7 


green solid circle 






2.0 - 6.3 


5 LMC regions 

Parker et al. (2001) 






+2.3 + 0.2 


blue solid triangle 






5-60 


NGC 1818 in LMC 








Santiago et al. (2001). outer region 




+2.5 




blue solid triangle 
NGC 1805 in LMC 




0.9-3 




Santiago et al. (2001). outer region 
blue solid triangle 




+3.4 
0.9-3 





Table 1. continued 
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a 


a 




r 71 jr l 

mass range [M & \ 


r 71 /r 1 

mass range \Mq\ 


mass range [M©J 


30 Dor* in LMC 








Selman et al. (1999). r > 3.6 oc 






+2.37 + 0.08 


cyan small open triangle 






3-120 


Selman et al. (19991. 1.1 < t/dc < 4.5 






+2.17 + 0.05 


cyan small open triangle 






2.8-120 


Sirianni et al. (20001 




+1.27+0.08 


+2.28 + 0.05 


cyan large open triangle 




1.35 - 2.1 


2.1 - 6.5 


Arches cluster* 








Figer et al. (19991. all radii 






+ 1.6 + 0.1 


cyan large solid circle 






6.3-125 


NGC 3603* 








Eisenhauer et al. (19981 




+1.73 


+2.7 


cyan small solid circle 




1 - 30 


15-70 


Globular clusters 








Piotto & Zoccali (19991 


+0.88 ±0.35 


+2.3 




yellow open triangles 


0.1 - 0.6 


0.6-0.8 




Galactic bulge 








Holtzman et al. ( 19981 


+0.9 


+2.2 




magenta filled square 


0.3 - 0.7 


0.7-1.0 




Zoccali et al. (20001 


+1.43 + 0.13 


+2.0 + 0.23 




magenta filled square 


0.15-0.5 


0.5-1.0 





Solar Neighbourhood (magenta dotted lines) 
Reid et al. (19991 +1.5 + 0.5 

0.02-0.08 

Herbst et al. (19991 < +0.8 

0.02 - 0.08 

Chabrier (2001. 20021 < +1 +1 / +2 
0.01 - 0.08 0.10 - 0.35 / 0.35 - 1.0 



Table 1. a(< Im >) data o btained since 1998 and until 2002. The data are shown in Fig. [201 in addition to the previously 
available data set compiled by Scalo (1998). Each a value is obtained at <lm>= (lm,2 — lm\)/2, Im = log 10 m, by the respective 
authors by fitting a power-law MF over the logarithmic mass range given by m\ and mi listed above. Some authors do not quote 
uncertainties on their a values. Notes: * are star-burst clusters; 1 thin green open circle emphasises the low-mass M35 datum 
which may be incomplete; 2 the mass range 1.35 < jti/Mq < 2.1 may be incomplete and is emphasised by the cross through the 
cyan large open triangle. 
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mass range 




r)N 






Vm 




P st 




[M Q ] 


[per cent] 






[per cent] 




[M /pc 3 ] 


[M /pc 2 ] 






Ct 3 






a 3 




a 3 


a 3 




2.3 


2.7 


4.5 


2.3 


2.7 


4.5 


4.5 


4.5 


0.01-0.08 


37.15 


37.69 


38.63 


4.08 


5.39 


7.39 


3.21 x 10" 3 


1.60 


0.08-0.5 


47.81 


48.50 


49.71 


26.61 


35.16 


48.21 


2.09 x 10~ 2 


10.45 


0.5-1 


8.94 


9.07 


9.30 


16.13 


21.31 


29.22 


1.27 x 10~ 2 


6.35 


1 - 8 


5.70 


4.60 


2.36 


32.38 


30.30 


15.09 


6.54 x 10~ 3 


1.18 


q 1 on 
o — 1ZU 


0.40 


0.14 


U.UU 


20.80 


7.83 


0.08 


3.63 x 10" 5 


6.53 x 10" 3 


to/M© = 


0.380 


0.292 


0.218 




P&t = 0.043 


Kot = 19.6 








2.3 


Oi 3 


= 2.7 




AAf c i/Af c i 




^max 


N cl 


M Gl 


N cl 


M cl 


m to 


[per 


cent] 








[Me] 




[M Q ] 


[M Q ] 


a 3 = 2.3 


a 3 = 2.7 




1 


16 


2.9 


21 


3.8 


80 


3.2 


0.7 




8 


245 


74 


725 


195 


60 


4.9 


1.1 




20 


806 


269 


3442 


967 


40 


7.5 


1.8 




40 


1984 


703 


1.1 x 10 4 


2302 


20 


13 


4.7 




60 


3361 


1225 


2.2 x 10 4 


6428 


8 


22 


8.0 




80 


4885 


1812 


3.6 x 10 4 


1.1 x 10 4 


3 


32 


15 




100 


6528 


2451 


5.3 x 10 4 


1.5 x 10 4 


1 


44 


29 




120 


8274 


3136 


7.2 x 10 4 


2.1 x 10 4 


0.7 


47 


33 



Table 2. The number fraction ?/jv = 100 J m ^ £(m) dm/ J m " £ (to) dm, and 
the mass fraction T]m = 100 f™*ml;(m)dm/M c i, M c \ = j™ u m£(m) dm, 
in per cent of BDs or main-sequence stars in the mass interval m\ to m^, 
and the stellar contribution, p st , to the Oort limit and to the Galactic-disk 
surface mass-density, E st = 2hp st , near to the Sun, taking mi = 0.01 Mq, 
m u = 120 Mf?) and the Galactic-disk scale-height h = 250 p c (m < 1 Mq 
iKroupa et all 1 19931) and h = 90 pc (to > 1 Mq, IScalol I1986D . Results are 
shown for the standard IMF (eq. |20|) . for the high-mass-star IMF approx- 
imately corrected for unresolved companions (a 3 = 2.7, to > 1M ), and 
for the PDMF (or 3 = 4.5, IScalol Il986t IKroupa et aHll993T ) which describes 
the distribution of stellar masses now populating the Galacti c disk. For 
gas S gas = 13 ± 3M Q /pc 2 and remnants S rom w 3M Q /pc 2 (|Weidemannl 
[1993). The average stellar mass is to = j™ u mti{m)dm/ f™ u l;(m)dm. 
N c i is the number of stars that have to form in a star cluster such that 
the most massive star in the population has the mass m max . The mass 
of this population is M c \, and the condition is J m £ i {m)dm — 1 with 

Jo ox™ £( w ) ^ TO = -^ cl — 1- ^M c i/M c i is the fraction of mass lost from the 
cluster due to stellar evolution, assuming that for to > 8 M Q all neutron stars 
and black holes are kicked out due to an asymmetrical sup ernova explosion, 
but that white dwarfs are retained (fWcidcma nn et al.lll992h and have masses 
towd = 0.084 mini + 0.444 [Mq]. This is a linear fit to the data in (jWeidemannl 
2000, their table 3) for progenitor masses 1 < m/M Q < 7 and towd = 0.5 Mq 
for 0.7 < to/Mq < 1. The evolution time for a star of mass mto to reach the 
turn-off age is available in Fig. [20j 





dN — £(to) dm — ^(m)dlm 
Cl(w) = (mlnl0)£(m) 

r M = *L (log 10 €L(im)) 
r = —x = 1 + 7 = 1 — a 

e.g. for power-law form: 


£ L =im r = Am- X 
£ = A 1 to" = A' to 7 
A' = A/lnlO 


gen 
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SalDeter ri955^ 
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L = -1.35 (a = 2.35) 
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thick long- dash- dotted line 
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Larson (1998) 

£/im short-dashed line 
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Hollenbach et al. (2005) 
not plotted in Fig. \20\ 


^L(m) = km 1 (1 — exp[— (TO,/TO c h) 7+1 ]) 
7 = 0.4, T = 1.35, to c1i = 0.18 M Q 




Roll 



Table 3. Summary of different proposed analytical IMF forms (the modern power-law form, the standard IMF, is presented 
in eq. [20)) . Notation: Im = log 10 (m/Af©) = ln(m./M©)/lnl0; dN is the number of single stars in the mass interval to to to + dm 
and in the logarithmic-mass interval Im to Im + dlm. The mass-dependent IMF indices, T(to) (eq. Gam), are plotted in Fig. 1201 
using the line-types defined here. Eq. MS was derived by Miller&Scalo assuming a consta nt star-format ion rate and a Galactic 
disk age of 12 Gyr (the uncertainty of which is indicated in the lower panel of Fig. 120k). iLarsonl (|1998l ) does not fit his forms 
(eqs. La and Lb) to solar-neighbourhood star-count data but rather uses these to discuss general aspects of likely systematic 
IMF evolution; the to q in eq. La and Lb given here are approximate eye-ball fits to the standard IMF. 
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